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The nucleon elastic electromagnetic form factors are 
fundamental quantities needed for an understanding of 
the nucleon's electromagnetic structure. The Sachs elec- 
tric, Ge, and magnetic, Gm, form factors [1], defined in 
terms of linear combinations of the Dirac and Pauli form 
factors, are of particular physical interest, as their evo- 
lution with Q^, the square of the four- momentum trans- 
fer, is related to the spatial distribution of charge and 
current within the nucleon. As such, precise measure- 
ments of these form factors over a wide range of are 
needed for a quantitative understanding of the electro- 
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magnetic structure not only of the nucleon, but also of 
nuclei (e.g., [2-4]). Further, in the low-energy regime of 
the nucleon ground state, the underlying theory of the 
strong interaction, Quantum Chromodynamics (QCD), 
cannot be solved perturbatively, and a proper description 
of even the static properties of the nucleon, the lowest 
stable mass excitation of the QCD vacuum, in terms of 
the QCD quark and gluon degrees of freedom still stands 
as one of the outstanding challenges of hadronic physics. 
Indeed, one of the most stringent tests to which non- 
perturbative QCD (as formulated on the lattice or in a 
model of confinement) can be subjected is the require- 
ment that the theory reproduce experimental data on 
the nucleon form factors (e.g., [5-7]). 

Because of the lack of a free neutron target, the neu- 
tron form factors are known with less precision than are 
the proton form factors, and measurements have been 
restricted to smaller ranges of Q^. A precise measure- 
ment of the neutron electric form factor, Geu, has proven 
to be especially elusive as the neutron's integral charge 
is zero. Prior to the realization of experimental tech- 
niques utilizing polarization degrees of freedom, values 
for Geu were extracted from measurements of the un- 
polarized quasielastic ^H(e,e'n)^H cross section and the 
deuteron elastic structure function A{Q'^). Those results 
for Geti deduced from measurements of the quasielas- 
tic ^H(e,e'n)^H cross section provided little information 
on Geu, as all results were consistent with zero over all 
ranges of accessed, < < 4 (GeV/c)^ (e.g., [8]). 
Similarly, results for Geu deduced from measurements of 
A(Q'^), although establishing G_e„ > for < < 0.7 
(GeV/c)^, were plagued with large theoretical uncertain- 
ties (~ ±40%) related to the choice of an appropriate 
A'' TV-potential for the deuteron wavefunction (e.g., [9]). 

With the advent of high duty-factor polarized elec- 
tron beam facilities and state-of-the-art polarized nu- 
clear targets and recoil nucleon polarimeters, experi- 
mental efforts over the past 15 years have now yielded 
the first precise determinations of Geu- Our experi- 
ment [10] was designed to extract the neutron electric 
to magnetic form factor ratio, GEn/GMn, from measure- 
ments of the neutron's recoil polarization in quasielastic 
^H(e, e'n)-'^H kinematics at three values of 0.45, 1.13, 
and 1.45 (GeV/c)^. These results were published rapidly 
by Madey et al. [11]; here we provide a more detailed 
report of the experiment and analysis procedures. 

The remainder of this paper is organized as follows. We 
begin, in Section II, with a brief overview of the exper- 
imental techniques utilizing polarization degrees of free- 
dom that have been employed for measurements of the 
neutron form factors. We continue with an overview of 
our experiment in Section III, and then discuss our neu- 
tron polarimeter in Section IV. Details of the analysis 
procedure are discussed in Section V. Our final results 
are then presented in Section VI and compared with se- 
lected theoretical model calculations of the nucleon form 
factors. Finally, we conclude with a brief summary in 
Section VII. A more detailed account of the discussion 



that follows may be found in [12]. 



II. NEUTRON FORM FACTORS 

A. Electron kinematics 

We will use the following notation for the electron kine- 
matics: (i!^e:Pe) wiU dcnotc the four-momentum of the 
initial electron, {Egi , pe' ) will denote the four- momentum 
of the scattered electron, 9e' will denote the electron scat- 
tering angle, uj = E^ — Eei will denote the energy transfer, 
q = Pe — Pe' will denote the three-momentum trans- 
fer, and Q2 = q2 _ ^ AE^Ee' sin^(6'e'/2) will denote 
the sqiiarc of the spacclike four-momentum transfer in 
the high-energy limit of massless electrons. The electron 
scattering plane is defined by Pe and Pe' . 



B. Measurements via polarized electron beams and 
recoil nucleon polarimetry 

1. Elastic N{e,e'N) scattering 

The polarization of the recoil nucleon, P, in elas- 
tic polarized-clcctron, unpolarized-nucleon scattering is 
well-known to be of the form [13-16] 



da 



= ao(pW+/iPW), 



(1) 



where ao denotes the unpolarized cross section, P^"^ de- 
notes the helicity-independent recoil polarization, P^'') 
denotes the helicity-dependent recoil polarization, and 
h = ±1 denotes the electron helicity. The polarization 
is customarily projected onto a {t,n,£) unit vector basis, 
with the longitudinal component, i, along the recoil nu- 
cleon's momentum; the normal component, n, perpendic- 
ular to the electron scattering plane; and the transverse 
component, perpendicular to the ^-component in the 
scattering plane. In the one-photon exchange approxi- 
mation, P^"^ ~ 0, and P'^''^ is confined to the scattering 



plane (i.e., P^'*' = 0). The transverse, P^"' , and lon- 
gitudinal, P(^\ components are expressed in terms of 
kinematics and nucleon form factors as [13-16] 



Be 



Pi'^) = Pe 



-2GeGm A/r(l + t) tan ^ 
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T-h 2r(H-r) tan^ 


^2 



(2b) 

where Pg denotes the electron beam polarization, r = 
(5^/4to^, and m denotes the nucleon mass. 



3 



Access to both P/''^ cx GeGm and P/''^ 



secondary analyzing reaction in a polarimcter is highly 
advantageous, as the analyzing power of the polarimeter, 
denoted Ay, and Pg cancel in the p/'^'/P^^'*'' ratio, yield- 
ing a measurement of Ge/Gm that is relatively insen- 
sitive to systematic uncertainties associated with these 
quantities. For the case of the neutron form factor ra- 
tio, as suggested by Arnold, Carlson, and Gross [16] and 
first implemented experimentally by Ostrick et al. [17], a 
vertical dipole field located ahead of a polarimeter config- 
ured to measure an up-down scattering asymmetry sen- 
sitive to the projection of the recoil polarization on the t- 
axis permits access to both P^'^^ and Pj''^ . During trans- 
port through the magnetic field, the recoil polarization 
vector will precess through some spin precession angle x 
in the t-t plane, leading to a scattering asymmetry, 
which is sensitive to a mixing of P/"' andPf\ 

C(x) = A,(p('')cosx + Pf^sinx) 
= A^|p('')|sin(x + (5). 



P 



(/i)x2i 1/2 



In the above, \P'^^^\ = [{Pi j ^ \^ i , , 
define the phase-shift parameter 5 according to 



(3) 



and we 



tan 6 ■ 



Gr 



cos ■ 



G 



M 



T + sin^ ■ 



(4) 



2. Quasielastic ^ H {e, e' nY H scattering 

The above formalism is directly applicable to an ex- 
traction of the proton form factor ratio, Gep/Gmp, from 
measurements of the proton's recoil polarization in elas- 
tic ^H(e, e'p) scattering. An extraction of the neutron 
form factor ratio, GEn/GMn, from measurements of the 
neutron's recoil polarization in quasielastic ^H(e, e'n)^H 
scattering is, however, complicated by nuclear physics 
effects, such as final-state interactions (FSI), meson ex- 
change currents (MEC), isobar configurations (IC), and 
the structure of the deuteron. The pioneering study of 
the sensitivity of the quasielastic ^H(e, e'n)'^H reaction 
to the neutron form factors, reported by Arenhovel [18], 
revealed that for perfect quasifree emission of the neu- 
tron (i.e., neutron emission along the three- momentum 

transfer q), p/''^ is proportional to Ge„, but is relatively 
insensitive to FSI, MEC, IC, and the choice of the NN- 
potential for the deuteron wavefunction. A more detailed 
study of the ^H(e, e'n)^H reaction reported by Arenhovel, 
Leidemann, and Tomusiak [19] found that these results 
also apply to P/^^ Similar findings were subsequently 
reported by [20, 21]. 

In Appendix A, we present a detailed discussion of the 
formalism for the kinematics and recoil polarization ob- 
servables for the quasielastic ^H(e, e'n)^H reaction. In 



particular, we provide there a definition for O^™', the 
polar angle between the proton momentum and q in the 
recoiling neutron-proton center of mass frame (hereafter, 
n-p cm. frame), a variable to which we will refer fre- 
quently throughout this paper. (Perfect quasifree emis- 
sion of the neutron is defined by 0^™' = 180°.) We 
follow this, in Appendix B, with a discussion of the sensi- 
tivity of the recoil polarization components to FSI, MEC, 
IC, and the choice of the A^A^-potential for the deuteron 
wavefunction at and away from perfect quasifree emis- 
sion. 



C. Measurements via polarized electron beams and 
polarized teirgets 

1. Elastic N {e, e' N) scattering 

The cross section in the one-photon exchange approx- 
imation for elastic polarized-electron, polarized-nucleon 
scattering is well-known to be of the form [13-15, 22] 



do- 
dOp' 



ao[l + hAeN{9*,(l)*)]. 



(5) 



Here, 9* and (j)* denote, respectively, the polar and 
azimuthal angle between the target nucleon polar- 
ization vector and q, and AeN{d* ,<!)*) denotes the 
polarized-electron, polarized-nucleon beam-target asym- 
metry, which is a function of kinematics and the nucleon 
form factors. The sensitivity of ^eJV to the form factors 
is enhanced if the target polarization is oriented in the 
electron scattering plane either parallel or perpendicular 
to q; in the former (latter) case, the expression for Aejv 
is identical to that for — P^''"' (Pt''*'') and will be denoted 
A|| {A±). Similar to the recoil polarization technique, 
measurements of both A± and j4|| are desirable as the 
target polarization (analog to the analyzing power) and 
beam polarization cancel in the ratio, again yield- 

ing a measurement of Ge/Gm that is relatively free of 
systematic uncertainties. 



2. Quasielastic ^H{e,e'nyH and^He{e,e'n) scattering 

The above formalism is directly applicable to a mea- 
surement of Gep/Gmp via the elastic ^H(e, e'p) reaction, 
but an extraction of Geu/Gmu from either the quasielas- 
tic ^H(e, e'n)^H reaction or the quasielastic ^He(e, e'n) 
reaction is again complicated by nuclear physics effects. 
For the case of the ^H(e, e'n)^H reaction, Cheung and 
Woloshyn [23] were the first to show that the polarized- 
electron, vector-polarized-deuterium beam-target asym- 
metry, A^^, is sensitive to Geh- More complete calcu- 
lations of AY^ that accounted for nuclear physics effects 
were later reported by Tomusiak and Arenhovel [24] and 
others [19, 21, 25, 26]. These calculations demonstrated 



4 



TABLE I: Chronological summary of published data on the neutron form factors from experiments employing polarization 
degrees of freedom and a recent analysis combining data on the deuteron quadrupole form factor, Gq, with data on t2o and 

T20- 



Reference 


Facility 


Published 


Type 


[(GeV/c)2] 


Quantities 


Note(s) 


Jones- Woodward et al. [31] 


MIT-Bates 


1991 


=^116(6, e') 


0.16 


A± Geu 


a, b 


Thompson et al. [32] 


MIT-Bates 


1992 


=^He(e, e') 


0.2 


Aj_, — > Gbji 


a, b 


Eden et al. [33] 


MIT-Bates 


1994 


2H(e,e'n) 


0.255 


J^t — *^En 


c,d 


Gao et al. [35] 


MIT-Bates 


1994 


^Ueie, e') 


0.19 


^11 Gmu 


a, e 


Meyerhoff et al. [37] 


MAMI 


1994 


^He(e, e'n) 
^He(e, e'n) 


0.31 


Aj_, A|| Geu 


a, b 


Becker et al. [38] 


MAMI 


1999 


0.40 


A_L, An Geu 


b,f 


Ostrick et al. [17], Herberg et al. [40] 


MAMI 


1999 


2H(e,e'n) 


0.15, 0.34 




b,c 


Passchier et al. [41] 


NIKHEF 


1999 


2H(e, e'n) 


0.21 


A^d * Geu 


6,c 


Rohe et al. [42], Bermuth et al. [43] 


MAMI 


1999/2003 


3He(e, e'n) 


0.67 


A^, Geu 


S,h 


Xu et al. [46] 


JLab 


2000/2003 


=^110(6, e') 


0.1 - 0.6 


A|| Gmu 


a, i 


Schiavilla and Sick [30] 




2001 


analysis 


0.00 - 1.65 


Gq — + Gsn 


3 


Zhu et al. [48] 


JLab 


2001 


2H(e,e'n) 


0.495 


Aed ~^ Geu 


6,c 


Madey et al. [11], this paper 


JLab 


2003 


2H(e,e'n) 


0.45, 1.13, 1.45 


p(h) p(h) ^ 

! — ^En 


c, k 


Warren et al. [50] 


JLab 


2004 


2H(e, e'n) 


0.5, 1.0 


^ed ~^ Geu 


c,9 


Glazier et al. [51] 


MAMI 


2005 


2H(e;e'n) 


0.30, 0.59, 0.79 


p(/t) p(/t) ^ 





"'Uncorrected for nuclear physics effects (i.e., for FSI, MEC, or IC). 
'Used the dipole parametrization for Gmu- 

Applied corrections for FSI, MEC, and IC by averaging calculations of Arenhovel et al. [18, 19, 24-26] over the acceptance. 
''Used the value for Gmu at = 0.255 (GeV/c)^ as measured by Markowitz et al. [34]. 
"^Used the Galster parametrization [36] for Geu- 
■'^Corrections for FSI and MEC calculated by Golak et al. [39]. 
^Used values for Gmu taken from the parametrization of Kubon et al. [44]. 

''Estimated corrections for FSI by scaling calculations of Golak et al. [45] at = 0.37 (GeV/c)^ to = 0.67 (GeV/c)^. 
'Used values for Geu taken from the parametrization of Hohler et al. [47]. 

^Theoretical analysis of data on the deuteron quadrupole form factor, Gq, tensor moment, t20, and tensor analyzing power, T20. 

*Used values for Gmu taken from the parametrization of Kelly [49]. 

'Used values for Gmu talsen from the parametrization of Friedrich and Walcher [52] . 



that for quasifrcc neutron kinematics, is strongly sen- 
sitive to Geu, but is relatively insensitive to FSI, MEC, 
IC, and the choice of the ATAT-potential for the deuteron 
wavefunction. 

For the case of the ^He(e, e'n) reaction, Blankleider 
and Woloshyn [27] were the first to study the sensitivity 
of the inclusive ^}i(2{e,e') asymmetry to Geu and Gmu- 
More detailed studios of the inclusive asymmetry carried 
out by others [28, 29] suggested that a clean extraction 
of Geti and Gun from the inclusive asymmetry would 
be extremely difficult due to proton contamination of the 
inclusive asymmetry. Such difficulties are, however, miti- 
gated in a ^Hc(e, e'n) coincidence experiment; as further 
motivation, Laget [21] demonstrated that the exclusive 
^He(e, e'n) asymmetry is relatively insensitive to the ef- 
fects of FSI and MEC for > 0.3 (GeV/c)^. 

D. Analysis of the deuteron quadrupole form factor 

The unpolarized elastic electron-dcuteron cross section 
is generally expressed in terms of the elastic structure 
functions, A{Q'^) and B{Q^). These are, in turn, func- 
tions of the deuteron's charge, Gq, quadrupole, Gq, and 
magnetic, Gm, form factors. Gq and Gq are of partic- 



ular interest for an extraction of Geti as they are both 
proportional to {Gep + Geu)- 

An unambiguous extraction of Gc, Gq, and Gm 
from a Roscnbluth separation of A{Q'^) and B{Q'^) re- 
quires some third observable. The tensor moments, t2j 
{j = 0,1,2), extracted from recoil polarization mea- 
surements in elastic unpolarized-electron, unpolarized- 
deuteron scattering, and the tensor analyzing powers, T2j 
{j = 0, 1, 2), as measured in elastic unpolarized-electron, 
polarized-deuteron scattering, are of particular interest 
as they are functions of Gq, Gq, and Gm [16, 22]. In- 
deed, after Gc, Gq, and Gm have been separated from 
A{Q'^), B{Q'^), and the polarization-dependent observ- 
ables, a value for Geu can be extracted from either Gq 
or Gq; however, as was shown by Schiavilla and Sick [30], 
an extraction of Geu from data on Gq is particularly 
advantageous as the contributions of theoretical uncer- 
tainties associated with short-range two-body exchange 
operators to Gq are small. 



E. Summary of results 

In Table I, we have compiled a complete chronologi- 
cal summary of all published data on the neutron form 
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FIG. 1: Current status of results for Geu ([H, 30, 33, 38, 40, 

41, 43, 48, 50, 51] and this work). The Galster parametriza- 
tion [36] is shown as the solid curve. See Table I for the 
reaction types for the individual data points. 



factors from experiments employing polarization degrees 
of freedom and a recent analysis combining data on the 
deuteron quadrupole form factor with the polarization- 
dependent observables t2o and T20. The ciirrent status of 
these results for Geu is shown in Fig. 1. We have omit- 
ted the results of Jones- Woodward et al. [31], Thompson 
et al. [32], and Meyerhoff et al [37] from this plot as these 
results were not corrected for nuclear physics effects. It 
should be noted that the results of Herberg et al. [40] 
and Bermuth et al. [43] supersede those of Ostrick et al. 
[17] and Rohe et al [42], respectively, as the former set 
reported the final results (corrected for nuclear physics 
effects) for their respective experiments. 

The range of Gsn is much more limited than those 
of the other three nucleon electromagnetic form factors, 
with only two results, those of Madey et al. [11] and the 
analysis results of Schiavilla and Sick [30] , extending into 
the > 1 (GeV/c)^ region. The agreement between 
these modern data and the Galster parametrization [36] 
with its original fitted parameters can be judged only as 
fortuitous. 



TABLE II: Nominal (central) values of the quasielastic elec- 
tron and neutron kinematics and neutron spin precession an- 
gles for each setting in the experiment. The data from the 
central values of 1.136 and 1.169 (GeV/c)^ were combined 
in our final analysis. 



Q2 




E,, 




T„ 


Precession 


[(GcV/c)^] 


[GeV] 


[GcV] 




[MeV] 


Angles X 


0.447 


0.884 


0.643 


52.65° 


239 


±40° 


1.136 


2.326 


1.718 


30.93° 


606 


0°,±90° 


1.169 


2.415 


1.789 


30.15° 


624 


±40° 


1.474 


3.395 


2.606 


23.55° 


786 


0°,±40°,±90° 



C specifically for this experiment, was used to measure 
the up-down scattering asymmetry arising from the pro- 
jection of the recoil neutrons' polarization on an axis per- 
pendicular to their momentum and parallel to the floor 
of Hall C. A vertical dipole field located ahead of NPOL 
was used to precess the recoil neutrons' polarization vec- 
tors through some chosen spin precession angle in order 
to measure this up-down scattering asymmetry from dif- 
ferent projections of the recoil polarization vector on the 
polarimeter's sensitive axis. This vertical dipole field also 
served as a sweeping field for the background flux of recoil 
protons from the deuteron target. 

Data were taken at four central values of 0.447, 
1.136, 1.169, and 1.474 {GeV/cf with associated elec- 
tron beam energies of 0.884, 2.326, 2.415, and 3.395 
GeV, respectively. The nominal (central) values of the 
quasielastic electron and neutron kinematics and the neu- 
tron spin precession angles, x, for each of these central 
Q"^ points are summarized in Table H. We note that the 
data acquired at the separate central Q"^ values of 1.136 
and 1.169 (GeV/c)^ were combined in our final analysis. 
Beam polarizations of 70-80% at currents of 20-70 /xA 
were typical throughout the duration of the experiment. 
The central axis of the neutron polarimeter was fixed at a 
scattering angle of 46.0° relative to the incident electron 
beamline for the duration of the experiment. The scat- 
tering asymmetries measured in our polarimeter were on 
the order of a few percent. 



III. EXPERIMENT 



Polarized electron source 



A. Overview of experiment 

Our experiment [10], E93-038, was conducted in Hall 
C of the Jefferson Laboratory ( JLab) during a run period 
lasting from September 2000 to April 2001. Longitudi- 
nally polarized electrons extracted from the JLab elec- 
tron accelerator [53] scattered from a liquid deuterium 
target mounted on the Hall C beamline. The scattered 
electrons were detected and momentum analyzed by the 
Hall C High Momentum Spectrometer (HMS) in coinci- 
dence with the recoil neutrons. A stand-alone neutron 
polarimeter (NPOL) [54], designed and installed in Hall 



Polarized electrons were produced at the accelerator 
source via optical illumination of a strained GaAs photo- 
cathode (GaAs on GaAsP [55]) with circularly polarized 
laser light from a high-power (~ 500 mW) Ti-sapphire 
laser [55, 56]; the linearly polarized light from the laser 
was circularly polarized with a Pockels cell. The helicity 
of the circularly polarized light emerging from the Pock- 
els cell was flipped at a frequency of 30 Hz (by switching 
the polarity of the high voltage applied to the Pockels 
cell) according to a pseudorandom scheme in which the 
helicity of one 33.3 ms window was randomly chosen, and 
the helicity of the following 33.3 ms window required to 
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be that of the opposite helicity (i.e., a sequence of such 
"hehcity pairs" could have been + — , — + , — h, H — , etc.). 
A A/2 plate was intermittently placed in the optics path 
upstream of the Pockels cell. This A/2 plate reversed the 
helicity of the electron beam that would otherwise have 
been induced by the Pockels cell, thereby providing the 
means for important systematic checks of any possible 
helicity-correlated differences. 



C. Hall C beamline 

Beam of the desired energy was extracted from the ac- 
celerator and then transported along the Hall C arc (se- 
ries of steering/bending magnets) and beamline. A num- 
ber of superharps [57] were used to monitor the beam 
profile, and four beam position monitors (cavities with 
four antennas oriented at angles of ±45° relative to the 
horizontal and vertical directions) provided absolute de- 
terminations of the beam position. The beam current was 
monitored with two monitors (cylindrical wave guides 
with wire loop antennas coupling to resonant modes of 
the beam cavity, yielding signals proportional to the cur- 
rent). 



D. Beam polarization measurements 

The beam polarization was measured periodically with 
a M0ller polarimeter [58] located along the Hall C beam- 
line approximately 30 m upstream of the cryotarget. We 
measured the beam polarization approximately every one 
to two days during stable accelerator operations. Mea- 
surements were also typically conducted following the in- 
sertion or removal of the A/2 plate at the polarized source 
or other major accelerator changes. A statistical preci- 
sion of < 1% was typically achieved after ~15-20 minutes 
of data taking. Details of the results of our beam polar- 
ization measurements will be discussed later, where it 
will be seen that the details of the analysis are relatively 
insensitive to the exact values of the beam polarization. 
Instead, the beam polarization information was primar- 
ily used to assess systematic uncertainties associated with 
temporal fluctuations in the polarization. 

It should be noted that although our production scat- 
tering asymmetry data were taken with beam currents 
as high as 70 /xA, the M0ller polarimeter was only de- 
signed for currents up to ~ 8 jiA. (due to the heating 
and subsequent depolarization of the iron target foil); 
therefore, it was necessary to assume that our beam po- 
larization measurements conducted at currents of 1-2 /xA 
were valid for the higher beam currents of our production 
running. The validity of this assumption has been veri- 
fied for operations in Hall A at JLab where the results of 
beam polarization measurements conducted at low cur- 
rents (M0llcr polarimeter) and high currents (Compton 
polarimeter) were found to agree [59]. 



E. Scattering chamber and cryotargets 

The scattering chamber consisted of a vertically- 
standing cylindrical aluminum chamber vacuum coupled 
to the incoming beamline. Two exit windows (made 
of beryllium) faced the HMS and NPOL, while an exit 
port faced the downstream beamline leading to the beam 
dump. During our experiment, the scattering chamber 
housed only one target ladder divided into a cryogenic 
target section and a solid target section. The cryogenic 
target section consisted of three cryogenic target "loops" . 
Each of these loops consisted of 4-cm and 15-cm long alu- 
minum target "cans", heat exchangers (heat loads from 
the electron beam were typically several hundred Watts), 
high- and low-power heaters (used to maintain the cry- 
otargets at their specified temperatures and to correct 
for fluctuations in the beam current), and various sen- 
sors. Liquid deuterium and liquid hydrogen, maintained 
at (nominal) operating temperatures of 22 K and 19 K, 
respectively, circulated through two of these loops; the 
third loop was filled with gaseous helium. Solid (car- 
bon) targets and 4-cm and 15-cm long "dummy targets" , 
composed of two aluminum foils spaced 4 cm and 15 cm 
apart, were mounted on the solid target section of the 
target ladder. As discussed in more detail later, data 
were taken with the dummy targets in order to assess the 
level of contamination due to scattering from the target 
cell windows. The thicknesses of the liquid deuterium 
and liquid hydrogen target cell windows were on the or- 
der of 4-6 mils, while those of the dummy targets were 
much thicker and on the order of 36-37 mils. 

To mitigate the effects of local boiling, the beam was 
rastered over a 2 x 2 mm^ spot on the cryotargets us- 
ing a fast raster system [60] located ~ 21 m upstream of 
the cryotargets. Target conditions (e.g., temperatures, 
heater power levels, etc.) were monitored continuously 
throughout the duration of the experiment using the 
standard Hall C cryotarget control system. 

F. High Momentum Spectrometer 

Scattered electrons were detected in the HMS, a 
three-quadrupole, single-dipolc (QQQD) spectrometer 
(all magnets are superconducting) with a solid angle 
acceptance of 6 msr (defined by an octogonally-shaped 
flared collimator), a maximum central momentum of 7.5 
GeV/c, a ±18% momentum acceptance, and a ~ 27 m 
flight path from the target to the detector package. 

1. Magnets 

The three quadrupolc magnets and the dipole magnet 
are mounted on a common carriage that rotates on a rail 
system about the target. The quadrupoles are 1.50 T 
maximum 20-ton (first, Qi) and 1.56 T maximum 30-ton 
(second, Q2, and third, Q3) superconducting coils with 
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FIG. 2: Schematic diagram of the ordering of the HMS detec- 
tor package elements. Shown axe the two drift chambers (DCl 
and DC2j, the two x-y hodoscopes (SIX/SIY and S2X/S2Y), 
the gas Cerenkov counter, and the lead-glass calorimeter. 



magnetic lengths of 1.89 m and 2.10 m, respectively. Qi 
and Qs are used for focusing in the dispersive direction, 
while Q2 provides transverse focusing. The dipole is a 
1.66 T maximum 470-ton superconducting magnet with 
a magnetic length of 5.26 m, a bend angle of 25°, and a 
bend radius of 12.06 m. 

The magnets were operated in their standard point- 
to-point tune in both the dispersive and non-dispersive 
directions. For our central points of 0.447, 1.136, 
1.169, and 1.474 (GeV/c)^, the nominal field strengths 
of Qi were 0.11, 0.31, 0.32, and 0.46 T; those of Q2 were 
0.13, 0.37, 0.38, and 0.55 T; those of Q3 were 0.06, 0.17, 
0.18, and 0.26 T; and, finally, those of the dipole were 
0.18, 0.47, 0.49, and 0.71 T. 



2. Detector package 

The detector package is enclosed within a concrete 
shielding hut and includes two drift chambers, two sets 
of hodoscopes, a gas Cerenkov counter, and a lead-glass 
calorimeter. A schematic diagram depicting the ordering 
of the detector package elements is shown in Fig. 2. 

a. Drift chambers The two multiwire drift chambers 
[61], used for tracking, each consist of six wire planes: (1) 
the X and X' planes, which provide position informa- 
tion on the x-coordinate (dispersive direction); (2) the 
Y and Y' planes, which provide position information on 
the y-coordinatc (non-dispersive direction); and (3) the 
U and V planes, which are inclined at ±15° angles rel- 
ative to the orientation of the X and X' planes. As 
seen by incoming particles, the ordering of these planes 
is XYUVY' X' . The active area of each plane is 113 (x) 
X 52 (y) cm^ with an alternating sequence of anode wires 
(25 /xm gold-plated tungsten) and cathode wires (150 /zm 
gold-plated copper-beryllium) spaced ~ 1 cm apart. The 
individual wire planes are separated by 1.8 cm, and the 
two drift chambers are separated by 81.2 cm. The cham- 
bers were filled with equal mixtures (by weight) of argon 
and ethane and maintained at a pressure slightly above 
atmospheric pressure. The signals from the anodes were 
read out in groups of 16 by multi-hit time-to-digital con- 
verters (TDCs). The fast branch of the signals from the 
hodoscope TDCs (to be described shortly) defined the 



TDC start for the electron arm trigger, while the de- 
layed signals from the drift chamber TDCs formed the 
TDC stop. 

b. Hodoscopes The x- (y-) planes of the two ho- 
doscopes, denoted S1X/S2X' (S1Y/S2Y), consist of 16 
(10) 75.5-cm (120.5-cm) long Bicron BC404 plastic scin- 
tillator bars with a thickness of 1.0 cm and width of 8.0 
cm. UVT lucite light guides and Philips XP2282B pho- 
tomultiplier tubes (PMTs) are coupled to both ends of 
each scintillator bar. The SlX/SlY and S2X/S2Y planes 
are separated by ~ 2.2 m. The fast branch of the PMT 
signals was routed to leading-edge discriminators. The 
discriminated signals were then split, with one set of out- 
puts directed to logic delay modules, TDCs, and scalers, 
and the other set directed to a logic module. The overall 
logic signaling a hit in any one of the hodoscope planes 
required a signal above threshold in at least one of the 
16 (10) PMTs mounted on the x > (y > 0) side of the 
bars and at least one of the 16 (10) PMTs mounted on 
the opposite .x < (y < 0) side. The slow branch of the 
PMT signals was directed to analog-to-digital converters 
(ADCs). 

c. Cerenkov detector The Cerenkov detector is a 
cylindrical tank (165-cm length and 150-cm inner diam- 
eter) filled with Perfluorobutane (C4F10, index of refrac- 
tion n = 1.00143 at STP). The pressure and temperature 
in the tank were monitored on an (approximately) daily 
basis and were observed to be highly stable. Pressures 
were typically ~0. 401-0. 415 atm (indices of refraction 
~1. 00057-1. 00059), translating into energy thresholds of 
~ 21 MeV (- 5.6 GcV) for pions (electrons). The tank 
is viewed by two mirrors, located at the rear of the tank, 
which focus the resulting Cerenkov light into two Burle 
8854 PMTs. The signals from these PMTs were directed 
to ADCs. During this experiment, information from the 
Cerenkov detector was used only for electron-hadron dis- 
crimination and not for HMS trigger logic purposes. 

d. Lead-glass calorimeter The calorimeter consists 
of 52 TFl lead-glass blocks stacked into four vertical 
layers of 13 blocks each. Each block has dimensions of 
70 X 10 X 10 cm^, corresponding to ~16 radiation lengths 
for the total four-laycr-tliic;kncss of 40 cm. As is indi- 
cated in Fig. 2, the four layers of the calorimeter are 
tilted at an angle of 5° relative to the central axis of the 
detector package to eliminate losses in the gaps between 
the individual blocks. Philips XP3462B PMTs are cou- 
pled to one end of each block, and the signals from these 
PMTs were routed to ADCs. Again, information from 
the lead-glass calorimeter was not used for HMS trigger 
logic purposes during this experiment. 



IV. NEUTRON POLARIMETER 

A. Overview 

A schematic diagram of the experimental arrangement 
with an isometric view of the neutron polarimeter is 
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FIG. 3: Isometric view of the NPOL flight path showing 
the Charybdis dipole magnet, the lead curtain, the front 

veto/taggcr array, the front array, the rear veto/tagger ar- 
ray, and the top and bottom rear arrays. 



shown in Fig. 3. The first element in the NPOL flight 
path was a dipole magnet (Charybdis) with a vertically 
oriented field that was used to precess the neutrons' spins 
through an angle x a horizontal plane. As a by- 
product, protons and other charged particles were swept 
from the acceptancx; during asymmetry measurements 
conducted with the field energized. The next item in the 
flight path was a 10.16-cm thick lead curtain, located di- 
rectly in front of a steel collimator (not shown in this 
figure). The lead curtain served to attenuate the flux of 
electromagnetic radiation and to degrade in energy the 
flux of charged particles incident on the polarimeter's de- 
tectors. 

The polarimeter consisted of 70 plastic scintillation 

detectors enclosed within a steel and concrete shielding 
hut. The front array of the polarimeter functioned as 
the polarization analyzer (via spin-dependent scattering 
from impolarized protons in hydrogen and carbon nuclei), 
while the top and bottom rear arrays, shielded by the col- 
limator from a direct line-of-sight to the target, were con- 
flgured for sensitivity to an up-down scattering asymme- 
try proportional to the projection of the recoil polariza- 
tion on a horizontally-oriented "sideways" axis (see next 
subsection). Double layers of thin- width "veto/tagger" 
detectors located directly ahead of and behind the front 
array tagged incoming and scattered charged particles. 
The flight path from the center of the target to the cen- 
ter of the front array was 7.0 m, and the distance from 
the center of the front array to the center of the rear 
array (along the polarimeter's central axis) was ~ 2.5 m. 



B. Polarimetry 

1. Coordinate systems 

Here we establish some necessary notation for a num- 
ber of diff'erent coordinate systems to which we will refer 
throughout the remainder of this paper. 

First, calculations of recoil polarization for the 
quasielastic ^H(e, e'n)^H reaction are usually referred to 
a {t, h, C.) reaction basis, defined on an event-by-event ba- 
sis in the n-p cm. frame according to 

I Hp--, n II q-- X p--, i = nx£, (6) 

where p^"^' and cf ™- denote, respectively, the incident 
neutron's momentum and the momentum transfer in the 
n-p cm. frame. The reaction basis can best be visualized 
by referring to the schematic diagram of the kinematics 
in the n-p cm. frame shown in Fig. 30 of Appendix A. 

Second, we define a polarim,eter basis, 
(^npol,^npol,5npol), fixed for all events, defined 
in the laboratory frame according to 

^NPOL II NPOL central axis, (7a) 
^NPOL ± Hall C floor, (7b) 

^NPOL = yNPOL X -2NPOL, (7c) 

with the center of the target defined to be the origin of 
this coordinate system. 

Third, the symmetric geometric configuration of the 
polarimeter's top/bottom rear arrays suggests the in- 
troduction of a polarimeter momentum basis, {S,N,L), 
which we again define on an event-by- event basis in the 
laboratory frame according to 

L II Pn, S II yNPOL XPn, N = L X S, (8) 

where p„ denotes a unit vector along the incident neu- 
tron's momentum in the laboratory frame. We will 
henceforth refer to the 5* and L axes as the polarimeter's 
"sideways" and "longitudinal'' axes of sensitivity, respec- 
tively. We express the recoil polarization in terms of the 
polarimeter momentum basis as P = PsS + PjyN + PlL. 

A schematic diagram showing the orientation of the 
polarimeter basis and polarimeter momentum basis co- 
ordinate systems is shown in Fig. 4. 

2. Detected scattering asymmetry 

We define NPOL polar and azimuthal scattering an- 
gles, denoted ^scat and ^scati according to 

sin^scat = \PnXp'n\, (9a) 
cos ^!>scat = S -u, (9b) 

where p'^ is a unit vector along the scattered neutron's 
three-momentum, and the unit vector u is defined ac- 
cording tOU= {pn X p'^)/\pn X p'J. 
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FIG. 4: Schematic diagram of the (4npol, ^npol, znpol) po- 
larimctcr basis (fixed for all events) and the (S, N, L) po- 
larimeter momentum basis (defined on an event-by-event ba- 
sis). Note that as pn is not, in general, restricted to the 
yNPOL-^NFOL plane, S is not, in general, parallel to xnpol- 



The cross section for elastic polarized-nucleon, 
unpolarized-nucleon scattering, denoted cr(0scat, "/"scat) 
for short, is of the form [62] 

(7(^scat,0scat) = Cro(6'scat) [1 + Ay{9sca.t)T' ■ u] 

« 0-o(6'scat) [1 + (6'scat )-Ps COS 0scat] , 

(10) 

where <Jo(0scat) and Ay [6 scat) denote the unpolarizcd 
cross section and the analyzing power, respectively. The 
above approximation is valid in the limit that Pm is small. 
It is then clear that the asymmetry, ^(^scat, '?i'scat), be- 
tween scattering "up" (S' • it < => cos (/)scat < 0) and 
scattering "down" (5 • it > =5> cos 0scat > 0) into in- 
finitesimal solid angles (6'scat, (/-scat) and (6'scat, <?i'scat +7r), 
respectively, for a particular value of Pg is 

i:(ft A> — 0'(^scat, '/'scat) ^ Cr(^scat, '/'scat + ■ 

sl.t'scat, '/'scat) — 7 \ '. 771 ) i 

O'l.fscat, 'Pscat) + Cr (^C/gcat, 

= (6'scat )^'S cos (/)sca 



/'scat + Tt) 



cat ■ 



(11) 



A single value of Ps is not, of course, presented to the 
polarimeter. Also, the top and and bottom rear arrays 
have a finite geometry; therefore, if the polarimeter is ge- 
ometrically symmetric in i^scat (i-C, geometrically sym- 
metric top and bottom rear arrays), the detected scat- 
tering asymmetry (i.e., averaged over kinematics and the 
top/bottom finite geometry), (^), is 



(0 - {Ps)Al' 



cff 



(12) 



where (Ps) and Ay^ denote, respectively, the acceptance- 
averaged value of the sideways component of the polariza- 
tion and the polarimeter's effective analyzing power av- 
eraged over its geometric acceptance (i.e., over cos(/)scat)- 
Henceforth, when we refer to the analyzing power Ay, it 
should be understood that we are referring to Ay^ . 



C. Charybdis dipole magnet and spin precession 

The Charybdis magnet was a water-cooled, 38-ton, 1.5- 
m tall, 2.3-m wide, and 1.7-m long iron dipole magnet 



TABLE III: Summary of the nominal values of the field in- 
tegrals (along the central axis) for the spin precession angles 
at each setting. P„ denotes the neutron velocity for the 
nominal (central) kinematics. 



Central 
[{GcY/cf] 


/3n 


Precession 
Angle X 


/|B| di 

[T-m] 


0.447 


0.604 


±40° 


0.6884 


1.136 


0.794 


±90° 


2.0394 


1.169 


0.799 


±40° 


0.9123 


1.474 


0.839 


±40° 


0.9576 


1.474 


0.839 


±90° 


2.1547 



installed in Hall C specifically for this experiment. The 
magnet was configured such that the gap between the 
pole pieces was 8.25 inches, and the geometric center 
of the magnet was located a distance of 2.107 m from 
the center of the target. The two poles were wired in 
parallel and powered with a 160 V-1000 A power sup- 
ply. Two-inch thick iron field clamps with apertures ma- 
chined to match the 8.25-inch pole gap were placed at 
the entrance and exit apertures resulting in an effective 
magnetic length of ~ 1.7 m. 

Calculations of the Charybdis field profile were per- 
formed with the TOSCA program [63] for various currents, 
and values for the field integral, / [B] dl, along the cen- 
tral axis were derived from these calctilations. The cur- 
rents were tuned for the various spin precession angles, 
X, according to the relation 



X 



Pn 



I 



B| d£. 



(13) 



where /Ujv is the nuclear magneton, g/2 — —1.913 for 
the neutron, and /3„ denotes the neutron's velocity. The 
field integrals for the precession angles at each of our 
points are tabulated in Table HI. 

The field along the central axis was mapped [64] at the 
conclusion of the experiment. We found that the values 
for the field integrals derived from our mapping results 
and the TOSCA calculations agreed to better than 0.76% 
for X = ±40° precession at = 0.447 {GeV/cf, 0.21% 
for X = +40° precession at = 1.169 (GeV/c)^, and 
0.35% for X = +40° precession at = 1.474 (GeV/c)2. 
Small differences in the measured field integrals for the 
two magnet polarities (corresponding to a ±0.3° spread) 
were observed for x = ±40° precession at = 0.447 
(GeV/c)^. Although we did not conduct field measure- 
ments for both polarities at the other points, it is 
reasonable to assume that the magnet behaved similarly 
for other current settings. 



D. Neutron polsirimeter physical acceptance 

The physical acceptance of the polarimeter was defined 

by a steel collimator with entrance and exit apertures 
located 483.92 cm and 616.00 cm, respectively, from the 
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FIG. 5: Schematic diagram (side view) of the NPOL shielding hut. The physical acceptance of the polarimeter, as defined by 
the collimator, is indicated by the dashed lines originating in the target. The rear array detectors were shielded from a direct 
line-of-sight to the target. The shadow shield, when inserted, was used to assess the room background rates. 



center of the target. The collimator was tapered, with the 
entrance (exit) port spanning a width of 72.6 cm (92.4 
cm) and a height of 37.3 cm (47.5 cm). The 10.16-cm 
thick lead curtain was located immediately upstream of 
the collimator's entrance port. 

A schematic diagram of the polarimctcr's shielding hut 
showing the shielding of the rear array detectors by the 
collimator from a direct line-of-sight to the target appears 
in Fig. 5. 

E. Neutron polarimeter detectors 

The polarimeter consisted of a total of 70 mean-timed 
BICRON-400 plastic scintillation detectors subdivided 
into a front veto/tagger array, a front array, a rear 
veto/tagger array, and symmetric top and bottom rear 
arrays. The front wall of the polarimeter's shielding hut 
was composed of 132.08-cm thick steel blocks; the only 
opening in this wall was the lead-shielded collimator. A 
schematic diagram of the polarimeter's detector configu- 
ration is shown in Fig. 6. 

1. Front veto/tagger array 

The function of the first series of detectors in the neu- 
tron flight path, the front veto/tagger array, was to iden- 
tify charged particles incident on the polarimeter. This 
veto array consisted of two vertically-stacked layers of five 
160.0 X 11.0 X 0.635 cm'' scintillators stacked with their 
long (160.0 cm) axes oriented horizontally and perpen- 
dicular to the central flight path and the thin (0.635 cm) 
dimension oriented along the flight path. The vertical 
spacing between the detectors in each layer was ~ 1 mm; 



therefore, to eliminate charged particle leakage, the two 
layers were offset from each other in the vertical direc- 
tion by ^ 1 cm. Each scintillator bar was coupled to two 
Philips XP2262 2-inch PMTs via plexiglass light guides. 



2. Front array 

The front array was segmented into 20 100 x 10 x 10 

cm^ scintillators; segmentation of the front array permit- 
ted us to run with luminosities as high as 3 x 10"'® cm~^ 
s~^ (70 /iA cmTcnt on a 15-cm liquid deuterium target). 
The long (100 cm) axes of these detectors were oriented 
horizontally and perpendicular to the central flight path 
and were stacked vertically into four layers of five de- 
tectors. The long ends of each scintillator were coupled 
via plexiglass light guides to 2-inch Hamamatsu R1828- 
01 PMTs powered by bases designed specifically for this 
experiment for purposes of high gain and highly linear 
output under conditions of high rate [65]. 



3. Rear veto/tagger array 

Similar to the front veto/tagger array, the purpose of 
the rear veto/tagger array was to identify charged parti- 
cles (e.g., recoil protons from np interactions in the front 
array) exiting the front array. The detectors in this ar- 
ray were identical to those in the front veto/tagger array 
and were vertically stacked in a similar fashion into two 
layers of eight detectors each. [We note that only one 
layer of eight detectors existed for the early part of the 
experiment during our = 1.136 (GeV/c)^ run).] As in 
the front veto/tagger array, each scintillator was coupled 
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FIG. 6: Schematic diagram (side view) of the NPOL detector configuration showing the top and bottom rear subarrays for 
measurement of an up-down scattering asymmetry. 



to two 2-inch Philips XP2262 PMTs. 



4- Rear array 

The top and bottom rear arrays each consisted of 
twelve detectors stacked into three layers of four detec- 
tors each. Each layer contained two "10-inch" 25.4 x 
10.16 X 101.6 cm"^ detectors sandwiched in between two 
larger "20-inch" 50.8 x 10.16 x 101.6 cm^ detectors. These 
detectors were oriented with their long (101.6 cm) axes 
parallel to the central flight path and their 50.8 cm or 
25.4 cm dimensions oriented horizontally. The centers of 
the inner, middle, and outer layers were located a verti- 
cal distance of 56.72 cm, 73.23 cm, and 89.74 cm, respec- 
tively, above or below the central axis of the polarimctcr 
and a horizontal distance of 2.52 m, 2.57 m, and 2.52 m, 
respectively, from the front array geometric center (see 
Fig. 6). The long ends of each scintillator were coupled 
via plexiglass light guides to 5-inch Hamamatsu R1250 
PMTs powered by the same bases built for the front ar- 
ray. 

The vertical positions of the top and bottom arrays rel- 
ative to the polarimeter's central axis were optimized for 
front-to-rcar scattering angles near the peak of the ana- 
lyzing power for np scattering (~ 15°-20° for our range 
of neutron energies). This configuration with scattering 
angles in the vicinity of ~ 15°-20° also guaranteed, for 
our kinematics, that only one of the nucleons (for elastic 
np interactions in the front array and assuming straight- 
line trajectories for the recoil proton through the front 



array) scattered into either the top or bottom array. We 
also note that the horizontal position of the middle de- 
tector plane was staggered relative to those of the inner 
and outer layers so that the majority of the front-to-rcar 
tracks passed through at least two of the three horizontal 
planes, reducing the dependence of the rear array detec- 
tion eflaciency on the scattering angle. 



F. Electronics, event logic, and data acquisition 

1. Electronics 

The signals from the 140 NPOL PMTs were processed 
with electronics sited in two locations: (1) one set, lo- 
cated inside the shielding hut was used to form the timing 
logic signal for each PMT (past experience with neutron 
time-of-flight and polarimetry experiments [66] revealed 
that locating the discriminators as close to the PMTs 
as practical yielded the best timing resolution); and (2) 
another set, located in the counting house, was used to 
define the logic for the various event types. 

A schematic diagram of the configuration of the elec- 
tronics in the shielding hut for each scintillator bar in the 
front and rear arrays is shown in Fig. 7. High voltage was 
applied to each PMT remotely by an EPICS-controlled 
64-channel high-voltage CAEN mainframe crate located 
in the counting house. Modest levels of high voltage were 
applied to the PMTs for the front array detectors, as 
deterioration in the performance of these PMTs was of 
concern because of the high count rates in these scintilla- 
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FIG. 7: Schematic diagram of the configuration of the electronics in the shielding hut for the front and rear array detectors. 
Note that the anode signals from the rear array detectors were not preamplified. 



tors; however, no deterioriation in their performance was 
observed during the experiment (instead, gains were sta- 
ble to within ~ 10%). To compensate for the resulting 
lower levels of gain obtained directly from these PMTs. 
the anode signals were preamplified by fast preamplifiers 
with a gain of eight, custom-designed and assembled for 
this experiment. The anode signals from the PMTs in 
the rear array and the front and rear veto/tagger arrays 
were not preamplified. 

The anode signals from the front and rear arrays were 
then directed to an LED driver and pulse height monitor. 
When desired, this device was used to assess the response 
of each PMT to a flashing bhic LED mounted on its light 
guide. The ccntroid channels of the LED spectra were 
monitored periodically, and any necessary changes to the 
high voltage levels were performed remotely. The gains 
of the front and rear veto/tagger array PMTs were not 
monitored with this system. 

The anode signals from all four detector arrays were 
then split. The signals in the fast branch (for the 
event trigger and timing measurements) were directed 
to either constant-fraction discriminators (front and rear 
arrays) or leading-edge discriminators (front and rear 
veto/taggcr arrays) located inside the shielding hut and 
then sent to the electronics in the counting house. We 
did not employ constant-fraction discrimination for the 
veto/taggcr array detectors for the following reasons: (1) 
the dynamic range of energy deposition in these detectors 



was small for those events of interest, so the time- walk 
was tolerable; and (2) the timing measurements from 
these detectors were not used for energy determinations, 
so rcsohitions of a few ns were sufficient for charged par- 
ticle tagging. Those signals diverted to the slow branch 
were routed through delays located inside the shielding 
hut and then sent to the counting house. 

Upon arrival in the counting house, both the ana- 
log and timing signals were directed through fil- 
ters/transformers designed to eliminate low- frequency 
noise. The analog signals were then sent directly to 
ADCs, while the timing signals were first sent to dis- 
criminators and then routed to two branches of a timing 
circuit. In one branch, the output from these discrimina- 
tors were directed through level translators, delays, dis- 
criminators, and then further split and directed to TDCs 
and scalers. In the other branch of this timing circuit 
(used to form the event triggers), the timing signals from 
the PMTs on all of the detectors, except those in the rear 
veto/tagger array, were first sent to logic modules which 
were used to generate logic signals for coincidences be- 
tween the timing signals for the two PMTs on each de- 
tector. Logical ORs were generated for each of the twenty 
front array detector two-PMT coincidences. These sig- 
nals were then sent to a fan-in with one set of outputs 
directed to scalers and the other through a discrimina- 
tor; the output from this discriminator was then directed 
to the trigger circuit. The logical ORs for the rear array 
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detectors and the front veto/tagger- array detectors were 
routed through a fan-in and then directed to the trigger 
circuit. The timing signals from the rear veto/tagger- 
array detectors were not used for trigger purposes. 



2. Event logic and triggers 

All event trigger logic was performed by two LeCroy 
8LM 2365 Octal Logic Matrix modules. Pretrigger logic 
signals from the HMS (coincident hits in at least three 
of the four hodoscope planes), the NPOL front array, 
the NPOL rear array, and the NPOL front veto/tagger 
array were routed to the 8LM modules. In addition to 
these logic signals, triggers from the polarized electron 
source were also input to these modules. As previously 
discussed, the helicity of the electron beam was flipped 
pseudorandomly at 30 Hz. Electronics at the polarized 
source generated a logic signal for readout of helicity- 
gated scalers for each 33.3 ms helicity window. Fur- 
ther, these modules also generated a helicity-transition 
logic signal which was used to veto otherwise valid data 
triggers that occured during transitions at the polarized 
source from one helicity state to another. The duration of 
this helicity-transition logic pulse was ~ 600 /xs, resulting 
in an effective data-taking helicity window of ~ 32.7 ms. 

An electronic module known as the Trigger Supervisor 
(TS) functioned as the interface between the 8LM logic 
modules and the data acquisition system (DAQ). The TS 
generated a logic signal indicating the status of the DAQ 
(e.g., busy or not busy) that was input to the logic mod- 
ules. The logic modules then determined whether the 
logic for any of the eight possible physics triggers (e.g., 
electron singles, electron/front array coincidences, elec- 
tron/front array/rear array coincidences, etc.) was satis- 
fied. If the logic for any particular trigger was satisfied, 
the TS generated an accept signal leading to generation 
of the appropriate ADC gate and TDC common signals. 
The ADCs, TDCs, and scalers were then read out with 
real-time UNIX-based processors. 

The event triggers of interest were three-fold coinci- 
dences between hits in the electron arm, the front array, 
and the rear array. These events constituted ~ 80-85% of 
the event triggers, as the higher rate events, such as elec- 
tron singles or two-fold coincidences between the electron 
arm and the front array, were prescaled. 



3. Data acquisition 

The DAQ was controlled by the CEBAF Onhne Data 
Acquisition System (CODA) [67]. CODA includes an 
event-builder subsystem programmed to assemble the in- 
dividual ADC channel, TDC channel, and scaler read-out 
data fragments into an event. The data for the events 
were then written to disk in CODA format by another 
subsystem. 



Typical data acquisition rates were one million events 
in ~ 1.0 (~ 0.5) hours with the Charybdis dipole field 
energized (de-energized). 



V. DATA ANALYSIS 
A. Electron reconstruction and tracking 

1. Overview of analysis code 

The raw ADC, TDC, and scaler data written to disk 
and encoded by the DAQ in CODA format were de- 
coded with a modified version of the standard Hall C 
ENGINE analysis code (see, e.g., [68] for a discussion of 
the standard version) employed for the analysis of nearly 
all experiments conducted in Hall C. Modifications to 
the standard version were necessary to accommodate the 
raw data stream from the 70 NPOL detectors; hereafter, 
whenever we refer to the ENGINE analysis code, it should 
be assumed that we are referring to our modified version 
of this code. 

For each event, the scattered electron's track through 
the HMS was reconstructed, and various kinematic quan- 
tities (e.g., momentum, energy, focal plane distributions, 
etc.) were computed. ENGINE was not configured to re- 
construct the track of the nucleon through the polarime- 
ter; instead, the NPOL detector data were simply writ- 
ten to new data files for later processing by other analysis 
tools. 



2. Extraction of electron information 

a. Tracking The overall strategy of the tracking al- 
gorithm [68] was to use the hit information from the drift 
chambers and reference start times provided by TDC in- 
formation from the scintillators in the hodoscope planes 
to reconstruct the tra)ec;tory of the; particle through the 
drift chambers. TDC information from those scintillators 
in the hodoscope planes recording hits was used to estab- 
lish reference start times. This information, coupled with 
TDC information from the drift chambers, was then used 
to determine the location of the hit in the drift chamber 
planes. "Left-right ambiguities" in the drift chambers 
(i.e., whether a particle passed to the left or right of any 
given wire) were resolved by fitting a (straight-line) track 
to each left-right hit combination in the six planes of each 
drift chamber. The full track through both drift cham- 
bers with the overall smallest track reconstruction was 
defined to be the final reconstructed track through the 
drift chamber planes. 

h. Transport ENGINE then attempted to relate the 
positions and angles at the focal plane (determined from 
the track through the drift chambers) to target quan- 
tities. In standard coordinate notation for transport 
through a spectrometer, Zf-p is taken to point along the 
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central ray of the spectrometer, ifp in the dispersive di- 
rection (by convention, taken to point "downwards"), 
and yip = zip x ifp. It should be noted that HMS focal 
plane variables are traditionally referred to the detector 
focal plane, defined to be perpendicular to the central 
ray (i.e., parallel to the drift chamber planes) with the 
origin of the xip-yip plane defined to be that point in 
space where the central ray of the spectrometer inter- 
sects the true (magnetic) focal plane. In addition to the 
dispersive and non-dispersive variables, two other stan- 
dard transport variables, and y^p, are defined to be 
the slopes of the rays at the focal plane, x'^^ = dxfp/dz 
and j/fp = dyip/dz, respectively. The focal plane vari- 
ables a:fp, yip, x^^, and y'f^ were converted to target 
quantities x[^^ = dxtar/dz, ytar, vLr = dytar/dz, and 
6 = (IPe'l — |p?"|)/|Pe^|, where |p^| denotes the central 
momentum setting, via computation of transport matrix 
elements derived from optics studies. For this choice of 
target coordinates, xtar was not reconstructed but was, 
instead, defined to be xtar — for all events. 



3. Sample electron reconstruction results 

Sample histograms of the reconstructed ^-distribution, 
hereafter referred to as the "Ap/p-distribution", at our 
lowest and highest points are shown in Fig. 8. The 
quasielastic peak is clearly visible in both spectra, but 
a large accompanying background of inelastic events as- 
sociated with pion-production in the target is present 
in the = 1.474 (GeV/c)^ spectrum. Inelastic peaks 
were also clearly visible in the — 1.136 and 1.169 
(GeV/c)^ spectra but are not shown here. A sample 
two-dimensional histogram of Ap/p plotted versus the in- 
variant mass, W, calculated from the electron kinematics 
according to 

W = V(^ + mA,)2-|q|2, (14) 



= 1 .474 {GeV/cf 




W [GeV/c^] 

FIG. 9: (Color online) Correlation plot of Ap/p versus W for 
the full HMS acceptance at = 1.474 (GeV/c)^ 



where mAr is the nucleon mass, is shown in Fig. 9 for our 
Q2 = 1.474 (GeV/c)2 point. The A(1232) resonance is 
prominent in this distribution. 

Hadrons in the HMS were identified via examination 
of the Cerenkov photoelectron spectrum. As expected, a 
hadron peak was not visible in the — 0.447 (GeV/c)^ 
spectrum; however, prominent hadron peaks (at zero 
photoelectrons) were observed at the three higher set- 
tings. An example of such a photoelectron spectrum from 
our Q2 = 1.474 (GeV/c)^ data is shown in Fig. 10. Guts 
on the number of photoelectrons, coupled with cuts on 
the energy deposition in the calorimeter, were sufficient 
for electron-hadron discrimination. 



B. Neutron polarimeter energy calibration 

The (charge-integrating) ADGs for the front and rear 
array detector PMTs were calibrated with the Gompton 
spectra from a "^^^Th source (2.61 MeV 7-rays); the front 
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FIG. 10: Cerenkov photoolcctron spectrum for the full HMS 
acceptance at = 1.474 (GeV/c)^. A prominent hadron 
peak appears at zero photoelectrons. 

and rear veto/tagger array detectors were not calibrated 
as ADC information was not used for charged particle 
tagging. These calibrations were parametrized in terms 
of an equivalent electron energy (denoted "eVee" ) , where 
the relation between the light output of recoil protons 
and Compton-scattered electrons in organic scintillator 
was found by Madey et al. [69] to be well described by 
the parametrization 

Te = ai [1 - cxp(-a2 [TpT] + a^Tp. (15) 

Here, Tp denotes the energy deposition of a recoil pro- 
ton, Te denotes the energy deposition of an electron that 
yields the equivalent light output, and the are empir- 
ically determined parameters. 

Unfortunately, the range of electron energies (2.38 
MeV Compton edge) was not sufficient, as typical energy 
depositions for the recoil protons were estimated to be > 
several MeVee [12]; further, the hardware thresholds for 
the front (rear) array detectors were set at 4 (10) McVce. 
To remedy these shortcomings, a custom-designed linear 
amplifier with a gain of ten was placed in the timing cir- 
cuit during calibration runs. The resulting ADC spectra 
were fitted to the sum of the Klcin-Nishina distribution 
(smeared by a Gaussian resolution function) and an ex- 
ponential background tail. Pulse- height calibrations were 
performed at three different times during the experiment 
(roughly at the start, middle, and conclusion); minor dif- 
ferences (~ 10%) in the extracted calibration parameters 
were observed but were deemed to be relatively unimpor- 
tant as the selection of quasielastic ^ll(e.e'n)^H events 
did not rely heavily on pulse height information. 



C. Neutron polarimeter timing calibration 

To optimize track reconstruction and background re- 
jection in the neutron polarimeter, the relative timing 



relationships between the NPOL detectors and the HMS 
were carefully calibrated with a scries of algorithms de- 
signed to: (1) generate position calibrations for each de- 
tector; (2) generate relative timing calibrations for each 
detector in the front array and discern the relationship 
between the mean time for each front array detector and 
the trigger mean time; (3) calibrate the timing between 
the HMS and the front array (yielding a coincidence time- 
of-flight); (4) generate relative timing calibrations for 
each detector in the rear array and calibrate the time- 
of-flight between the front array and the rear array; and 
(5) generate position and timing calibrations for the front 
and rear veto/tagger detectors. 



1. Front and rear array position calibrations 

The position calibration algorithm for the front and 
rear array detectors employed data acquired with the 
Charybdis magnet dc-encrgized, such that charged par- 
ticles illuminated the front array almost uniformly. The 
relationship between the hit position and the difference 
(in channels) between the TDCs from the PMTs mounted 
on the two ends of each scintillator was parametrized in 
a linear form with an unknown slope and offset. His- 
tograms of these TDC channel differences were accumu- 
lated for each detector and then boxcar-smoothed. The 
algorithm identified the channel of maximum content and 
then scanned away in both directions until channels with 
10% of the maximum content were identified. Slope and 
offset parameters were then chosen such that these 10%- 
content channels were aligned with the physical edges of 
each detector; the resulting calibrated position spectra 
displayed sharp edges near the physical detector edges. 



2. Front-array timing and trigger calibrations 

The first goal of the front-array timing calibration was 
to align the mean times of all the detectors in the front 
array using events with a single hit in the front array. 
Data acquired with the Charybdis magnet energized (for 
suppression of background processes) were employed for 
this step of the timing calibration, and events with > 
(> 1) hits in the front veto/tagger array (front array) 
were discarded. An offset was chosen for each detector 
such that the mean value of its mean-time spectrum was 
aligned on zero. 

The second goal of the front array timing calibration 
was to construct a variable that could be used to identify 
which hit generated the trigger (for events with multiple 
front array hits), as the trigger circuit did not identify 
the triggering hit. Proper identification of the triggering 
hit via examination of the correlation between the TDC 
channels for the two PMTs on each detector and the po- 
sition dependence of the mean times yielded self-timing 
spectra with FWHM of ~ 0.4 ns. 
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3. Coincidence time-of-flight calibrations 

To maximize our signal-to- noise ratio, we constructed 
a coincidence' timc-of-fliglit variable that accounted for 
the quasielastic ^H(e, e'n)^H kinematics, pathlength vari- 
ations through the HMS and NPOL, and variations in the 
delay between an interaction in a detector and the arrival 
of its timing signal at the TDC. For this step of the cali- 
bration, a minimal set of cuts were applied to the data for 
purposes of (loose) quasielastic event selection (e.g., cuts 
on the calorimeter energy deposition, Ap/p, etc.). Again, 
front array single-hit events (with no hits in the front 
veto/tagger array) acquired with the Charybdis magnet 
energized were used for this step of the calibration. 

The algorithm first predicted the neutron time-of-flight 
from the target to the front array using only position 
information (i.e., the reconstructed vertex information 
for the primary scattering event in the target cell and the 
position of the front array hit) and electron kinematics. 
For a three-body final state (i.e., no pion production), 
four-momentum conservation demands 

rud + oj = ViPnP + ™n + ^IPpP (16a) 
q = P™ + Pp. (16b) 

From this, it follows that a value for |p„| (and, then, the 
predicted neutron time-of-flight) can be derived from the 
solution to the quadratic equation A|p„|^ -|- B|p„| -|- C = 
0, where 

A = {md + oj)^ - {(l-pn)^, (17a) 
B = -2(q-p„)£), (17b) 
C = ml{md+u;f-D^, (17c) 
2D = ml + ml - ml - Q'^ + 2mdUJ. (17d) 

A value for the actual measured time-of-flight was then 
extracted from information in the signal output of a TDC 
started by a signal generated by the NPOL trigger and 
stopped by the HMS trigger, a correction for pathlength 
variations and delays between interactions and signals in 
the HMS computed by ENGINE, and the mean time of 
the front array detector recording the hit. This mea- 
sured time-of-flight was then compared with the pre- 
dicted time-of-flight, and the resulting difference, the co- 
incidence time-of-flight (hereafter, referred to as cTOF), 
was computed for each event. The resulting cTOF spec- 
tra were fairly narrow with FWHM of 1.25 ns and 
signal-to-noise ratios of ~ 6:1-10:1. Sample cTOF spec- 
tra are shown later in this paper. 

4- Rear-array timing calibrations 

The algorithm for the rear-array timing calibration se- 
lected single-hit events (with no hits in both the front 
and rear veto/tagger arrays) acquired with the Charyb- 
dis magnet energized and then filtered these hits accord- 
ing to a set of cuts designed to select quasielastic events. 
In addition, a |cTOF| < 2 ns cut was enforced. 



In the first step, the algorithm aligned the mean time 
spectra of the rear array detectors relative to each other. 
As for the front array, histograms of mean times were 
accumulated for each detector. The channel of maximum 
content was identified, and an offset parameter for each 
detector was then chosen such that the peak channel was 
aligned on zero. 

In the second step, the algorithm performed an abso- 
lute timing calibration of the rear array detectors relative 
to the front array detectors via a front-to-rear velocity 
calibration. The scattering angle for the front-to-rear 
track was computed using the incident neutron's three- 
momentum and the position information for the hits in 
the front and rear array. The algorithm then predicted 
the front-to-rear velocity for elastic np scattering in the 
front array via computation of the scattered neutron's 
kinetic energy, T„p, where 

rr. 2T„ cos^ ^scat , . 

(7„ + l)-(7n-l)cOs2escat" ^ ^ 

Here, T„ denotes the incident neutron's kinetic energy, 
^scat denotes the neutron scattering angle in the po- 
larimeter, 7„ is the usual Lorentz factor for the incident 
neutron, and the proton and neutron masses are assumed 
to be equal. Relative time-of-flight (hereafter, referred 
to as rTOF) histograms, deflned to be the difference be- 
tween the predicted and measured values of the front- 
to-rear time-of-flight, were accumulated, and offsets were 
then chosen for each detector such that the peak channel 
was aligned on zero. Again, sample rTOF spectra are 
shown later in this paper. 

5. Front and rear veto/tagger-array calibrations 

The position and timing calibration of the front 
and rear veto/tagger- array detectors consisted of three 
steps. Data for charged particle tracks acquired with 
the Charybdis magnet de-energized were employed for 
this calibration; hits were required in each layer of the 
front veto/tagger array, the front array, and the rear 
veto/tagger array. 

First, as leading-edge discrimination was employed for 
these detectors, the algorithm began by computing cor- 
rections for walk. The relationship between the ob- 
served TDC and ADC channels, TDCobs and ADCobs, was 
parametrized as TDCobs = TDC -|- 7log(ADCobs/ADCpeak), 
where TDC denotes the TDC channel in the absence of 
walk effects, 7 is an empirical parameter, and ADCpeak 
denotes the peak ADC channel. A value for 7 was then 
computed via the method of least squares. 

Second, the veto/tagger array detectors were position 
calibrated using a different algorithm than that employed 
for the position calibration of the front and rear array 
detectors due to the facts that the collimator partly ob- 
scured the edges of the front veto/tagger array detec- 
tors and that the outer rear veto/tagger array detectors 
did not receive adequate illumination from front-to-rear 
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charged tracks. (The front and rear veto/tagger arrays 
were designed to provide more than adequate coverage 
of target-to- front and front-to-rear charged tracks.) As 
such, position cahbration parameters for these detectors 
were deduced via a comparison of the recorded hit posi- 
tion with the nearest hit position in the front array, and 
offset parameters were determined via a minimization 
of the difference between the predicted and recorded hit 
positions. To improve the statistics for the outer rear 
array veto/tagger detectors, the algorithm searched for 
(n, p) charge-exchange events in the front array. Tracks 
from these events were used to predict hit locations in 
the rear veto/tagger array detectors, and position cali- 
bration parameters were then deduced from another 
minimization of the difference between the predicted and 
recorded hit positions. The resulting calibrated position 
spectra were well aligned about the physical center of 
each detector with somewhat more rounded spectra than 
observed in the front and rear array spectra due to the 
use of leading-edge discrimination. 

Last, the mean times were aligned relative to each 
other via the same procedure employed for the mean- 
time calibration of all the other detectors. 



D. Nucleon reconstruction and tracking 

1. Overview of analysis code 

The algorithm we developed for reconstruction and 
tracking in the neutron polarimeter began by translat- 
ing the raw NPOL detector data decoded by ENGINE into 
hit positions and times. The code then attempted to de- 
termine which hit in the front array generated the trigger. 
All hits were then filtered according to a number of dif- 
ferent selection criteria, with the surviving hits grouped 
into recognizable patterns. The code then attempted to 
determine the primary hits in the front and rear arrays 
and the charges of the incident particle and the particle 
detected in the rear array. Finally, kinematic quantities 
and timc-of-flight variables were then computed for those 
events satisfying all tracking criteria. 



2. Trigger selection and hit filtering 

The algorithm assigned the location of the triggering 
front array hit to the detector with the smallest absolute 
self timing value. All hits were then filtered according 
to a munber of selection criteria designed to discard hits 
with unphysical reconstructed detector positions or mean 
times falling outside of specified windows. These mean 
time windows were chosen sufficiently wide for purposes 
of quasielastic event selection, elastic/quasielastic scat- 
tering in the front array, and charged particle tagging 
in the veto/tagger arrays. In particular, the mean-time 
windows for both the front and rear veto/tagger arrays 



safely bracketed the entire peak regions with the borders 
extending into the regions of flat background. 



3. Pattern grouping and track reconstruction 

a. Incom,plete and simple events The algorithm be- 
gan by identifying incomplete and simple events. First, 
events with either no surviving hits in the front and/or 
rear array or events with hits in both the top and bottom 
rear array were discarded. Second, simple events with ex- 
actly one hit in the front array, one hit in the rear array, 
and no hits in both the front and rear veto/tagger arrays 
were identified. For these events, the incident particle 
and the particle detected in the rear array were, obvi- 
ously, designated neutral particles, and reconstruction of 
the track was deemed complete. 

b. Multiple hit events The majority of the events 
were more complicated than these simple events because 
of propagation of the recoil protons through adjacent 
scintillator bars or multiple scattering of the neutron. 
For these more complicated events, the code began by 
identifying which layer in the front array (i.e., first, sec- 
ond, third, or fourth) was hit first; henceforth, we will 
refer to the hit(s) in this layer as the "first cluster". If 
the first cluster contained more than one hit, the (ver- 
tically) highest and lowest hits were identified; such hit 
patterns were assumed to be the result of an np or pp 
interaction in one detector followed by the penetration 
of the recoil proton into a vertically adjacent detector. 
Accordingly, if the hits occurred in non-contiguous de- 
tectors within the same vertical layer (i.e., existence of a 
vertical "gap"), the event was discarded. 

The code then searched for evidence of one or more 
"missing layers" in the front array (e.g., an event with 
hits in the first layer and the fourth layer); a missing 
layer was taken to be evidence for multiple scattering of 
the incident neutron. If such a "second cluster" of hits 
was not found, the location of the front array scatter- 
ing vertex was assigned to the highest (lowest) hit in the 
first cluster if the top (bottom) rear array recorded one 
or more hits. If, instead, a second cluster of hits was 
foimd, the code determined whether the second cluster 
contained a gap; again, events with gaps in the second 
cluster were discarded. The algorithm then attempted 
to discern whether the second cluster was located above 
or below the first cluster; if the second cluster was above 
(below) the first cluster, the location of the first cluster 
scattering vertex was assigned to the highest (lowest) hit 
in the first cluster. Then, if the top (bottom) rear array 
was hit, the location of the second cluster scattering ver- 
tex was assigned to the highest (lowest) hit in the second 
cluster. Finally, if more than one hit was recorded in 
either the top or bottom rear array, the rear array scat- 
tering vertex was assigned to that hit closest in distance 
to the final front array scattering vertex. 

Illustrative examples of two possible types of recon- 
structed tracks are shown in Fig. 11. We note here, and 
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FIG. 11: Examples of reconstructed tracks for: (a) an event with a single cluster in the front array, no missing layers, and 
multiple hits in the top rear array; and (b) an event with two clusters in the front array (separated by two missing layers) and 
multiple hits in the bottom rear array. 



discuss later in Section VE2, that events with a "sec- 
ond cluster" were reconstructed but were not used in our 
extraction of scattering asymmetries. 

4- Charge identification 

After the track through the front and rear arrays was 
reconstructed, the code then checked for hits in the 
veto/tagger arrays. The charge of the incident parti- 
cle was determined via the following algorithm. (1) If 
there were no hits in any of the front veto/tagger detec- 
tors, the particle was designated a neutral particle. (2) 
If there were hits in the front veto/tagger detectors, the 
radial distance between the location of the veto/tagger 
hit and the location of the f irst scattering vertex was 
computed according to d = \/ (xvt — a^fr)^ + (j/vt — 
where the coordinates refer to the polarimctcr basis, de- 
fined in Eq. (7). If at least one hit in each veto/tagger 
layer satisfied d < 30 cm, the incident particle was desig- 
nated a charged particle. If no hits in cither veto/tagger 
layer satisfied d < 30 cm, the incident particle was des- 
ignated a neutral particle. Finally, if a hit in one of the 
front /veto tagger layers satisfied this distance require- 
ment but no hits in the other layer satisfied this condi- 
tion, the charge of the incident particle was declared to 
be ambiguous. 

The algorithm for the determination of the charge 
of the particle detected in the rear array was essen- 
tially identical to that described above. The only dif- 
ference was that the code predicted where the hits in the 
rear veto/tagger arrays should have occurred assuming a 
straight-line trajectory from the final front array scatter- 
ing vertex to the rear array scattering vertex. The com- 
puted value of the radial distance between the location of 
the actual hit and the predicted hit was then used, in an 
identical manner, for rear array neutral/charged tagging. 

The choice of the 30-cm radial track-distance threshold 
was based on an examination of track-distance spectra 
for the front and rear veto/tagger arrays. The spectra 
for the front veto/tagger array were found to be rela- 
tively narrow with an abrupt change in slope around 30 
cm, believed to be related to these scintillators' position 
resolution. The spectra for the rear veto/tagger array 



did not contain such a feature as the recoil protons aris- 
ing from interactions in the front array were widely dis- 
tributed in angle; nevertheless, the same 30-cm condition 
was employed as the position resolutions for these detec- 
tors were similar to those in the front veto/tagger array. 

5. Kinematic distributions and time-of-flight variables 

Following K^construction of the track through the po- 
larimctcr, kinematic and time-of-flight quantities were 
computed for fully reconstructed events. First, the in- 
cident particle's momentum was computed using only 
position information for the reconstructed target vertex, 
position information for the first scattering vertex in the 
front array, and the four- momentum transfer (w, q), via 
solution of the quadratic equation for |p„| given previ- 
ously in Eq. (17). The momentum was then used to pre- 
dict the target-to- front array time-of-flight; the difference 
between the predicted and measured time-of-flight was 
then stored as the cTOF variablc\ Laboratory frame po- 
lar and azimuthal neutron scattering angles with respect 
to q, 9nq and 4>nq, were computed from information on q 
and p„ . Second, front-to-rear polar and azimuthal scat- 
tering angles, (?scat and (/)scat, were computed using in- 
formation on p„ and the scattering vertices in the front 
and rear arrays. This information was used to compute 
a value for T„p, Eq. (18), which was then used to predict 
the front-to-rear time-of-flight; the difference between the 
predicted and measured time-of-flight was then stored 
as the rTOF variable. Finally, the missing momentum, 
Pmiss, missing energy, Ey^iss, and missing mass, rrimiss, 
were computed according to 

Pmiss = q - Pn, (19a) 

Emiss = {rud + Lo)- (T„ + m„), (19b) 

"Jmiss = \/i^miss- IPmissP- (19c) 

6. Sample nucleon reconstruction results 

To illustrate the full range of the polarimeter's accep- 
tance, sample two-dimensional histograms of |pmiss| plot- 
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FIG. 12: (Color online) Correlation plot of |pmiss| versus W for the full NPOL acceptance at = 1.136 and 1.474 (GeV/c) 



ted versus the invariant mass W at our = 1.136 and 
1.474 (GeV/c)^ points are shown in Fig. 12. A minimal 
set of cuts designed to eliminate scattering from the tar- 
get cell walls, hadrons in the HMS, and protons incident 
on NPOL were applied to these spectra. Our acceptance 
was sensitive to missing momenta ranging up to 450 
MeV/c at our highest point. As can clearly be seen 
in these correlation plots, quasielastic events were asso- 
ciated with missing momenta in the range < 150 MeV/c. 
Larger values of |pmiss| are, of course, seen to corre- 
spond to inelastic events, with the A(1232) resonance 
prominent at large missing momenta in the = 1.474 
(GeV/c)^ spectrum. The correlation plot for = 1.169 
(GeV/c)^ was essentially identical to that at — 1.136 
(GeV/c)2, while the = 0.447 (GeV/c)^ distribution 
was restricted to considerably smaller ranges of |pmiss| 
(< 100 MeV/c). 



E. Data selection criteria, data sets, and cuts 

1. Data selection criteria and data sets 



TABLE IV; Quantity of data (accumulated charge) employed 
for the final production analysis. A total of 194 Coulombs of 
charge was delivered to the experiment for production running 
with the deuterium target. 



Central 


Precession 


Charge 


KGeV/cf] 


Angle X 


[Coulombs] 


0.447 


-40° 


25.122 


0.447 


+40° 


14.569 


1.136 


0° 


27.587 


1.136 


-90° 


4.701 


1.136 


+90° 


4.158 


1.169 


-40° 


7.006 


1.169 


+40° 


6.321 


1.474 


0° 


26.239 


1.474 


-90° 


4.097 


1.474 


+90° 


4.098 


1.474 


-40° 


20.803 


1.474 


+40° 


16.762 


Total 




161.463 



neutron spin precession angles. 



Only those data runs satisfying the following criteria 
were employed for the final production data analysis: 
(1) no problems with the HMS equipment (e.g., mag- 
net trips, detector failures, etc.); (2) no problems with 
delivery of the electron beam (e.g., unstable beam pa- 
rameters); (3) no problems with the DAQ; (4) no prob- 
lems with the cryogenic target (e.g., large temperature 
fluctuations, monitoring system failures, etc.); and (5) 
no problems with the Charybdis magnet or the NPOL 
detectors (e.g., fluctuations in the magnet current, de- 
tector high-voltage trips, etc.). We note that additional 
problems may have resulted in the designation of a run 
as unsuitable for the production analysis. 

The quantity of data satisfying the above selection cri- 
teria is summarized in Table IV. There, we list the accu- 
mulated charge for each of the individual points and 



2. Cuts for extraction of time- of -flight spectra 

A summary of the final set of cuts applied to the pro- 
duction data sets for extraction of the cTOF and rTOF 
time-of-flight spectra is as follows. 

a. Target variables Scattering from the target cell 
windows was suppresed via the requirement that the re- 
constructed target vertex lie within ±7 cm of the center of 
the target (for the 15-cm target) along the incident beam- 
line. Further, events with unreasonable reconstructed 
values for xj^,. and y[^^ were discarded. 

h. HMS variables The reconstructed electron track 
was required to fall within the the collimator acceptance, 
and events with unreasonably large track reconstruction 
values were discarded. Hadrons in the HMS were 
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suppressed via cuts on the number of Cerenkov photo- 
electrons and the energy deposition in the calorimeter. 
Events away from the quasielastic peak were suppressed 
via a tight Ap/p G [-3%, +5%] cut. 

c. NPOL variables Software thresholds of 8 (20) 
MeVee designed to suppress low-energy backgrounds 
were applied to the front (rear) array pulse height dis- 
tributions. Also, to suppress lower-energy neutrons orig- 
inating from charge-exchange Pb(p, n) reactions in the 
lead curtain (discussed in more detail later), the mean 
times for front array hits were required to lie within a 
[—5,5] ns window, due to the expected degradation in 
the energy of the incident protons prior to the charge- 
exchange reaction. Events with more than one scatter- 
ing vertex in the front array (i.e., existence of a second 
cluster) were discarded to eliminate the effects of depo- 
larization following the first interaction in the front array. 

The front-to-rear polarimeter scattering angle, ^scatj 
was required to satisfy 6'scat G [5°, 35°] at = 0.447 
(GeV/c)2 and e [5°, 30°] for the other points. The 
lower cut of 5° eliminated unreasonably small scattering 
angles, while the upper cut of 30° or 35° was used to 
suppress zero (or negative) values of the analyzing power 
at larger scattering angles (as predicted by SAID [70]). 

d. ^H{e,e'n)^H reaction variables Pion-production 
events were suppressed via tight cuts on the missing mo- 
mentum and invariant mass of [pmissl < 100 MeV/c and 
W < 1.04 GeV/c2. 



F. Extraction of time-of-flight spectra and 
scattering asymmetries 

1. Polarimeter event types 

An analysis code developed to extract the physical 
scattering asymmetries subjected each event to the cuts 
discussed previously. In addition, each event was also 
subjected to a more stringent test for the determination 
of the incident particle's charge. As we used single-hit 
TDCs, an early accidental hit in a front veto/tagger de- 
tector falling outside the mean-time window for the front 
veto/tagger array would have prevented that TDC from 
recording any later (on-time) hits, leading to the incor- 
rect tagging of a charged particle as a neutral particle. 

Histograms of cTOF were accumulated for two types 
of front array scattering events, (n, n) and (n,p) events, 
corresponding (for a neutral particle incident on the po- 
larimeter) to the detection of a neutral and charged par- 
ticle, respectively, in the rear array. We identified (n, n) 
events with the scattering of the neutron from the front 
array to the rear array, while we identified {n,p) events 
with forward scattering of the recoil proton with suffi- 
cient energy for penetration of the front array. It should 
be noted that for the incident neutron kinetic energies 
of interest, the analyzing power for elastic np scattering 
becomes negative for neutron scattering angles greater 
than ~ 40°; therefore, the signs of the detected asymme- 
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FIG. 13: (Color online) Correlation between cTOF and iTOF 
at = 1.474 (GeV/c)^ with the various event types (see 
text) identified. 



tries for (n,n) and (n,p) events were the same. Events 
with charges deemed ambiguous in either the front or 
rear array were rejected. 

Histograms of rTOF summed over all front-to-rear 
tracks were accumulated for those events falling within a 
prescribed cTOF window. To compensate for variations 
in the flight path between the front array and the rear 
array, the rTOF values were normalized to a nominal 250- 
cm flight path. The accumulated rTOF spectra were de- 
composed into the following event types: (1) "RU events" 
(positive beam helicity and scattering from the front ar- 
ray to the top rear array); (2) "LU events" (negative 
beam helicity, top rear array); (3) "RD events" (positive 
beam helicity, bottom rear array); and (4) "LD events" 
(negative beam helicity, bottom rear array). The scat- 
tering asymmetries were then extracted from the yields 
in these four spectra. 



2. HMS-NPOL coincidence event types 

A two-dimensional histogram of the correlation be- 
tween cTOF and rTOF summed over (n, n) and {n^p) 
events at = 1.474 (GeV/c)^ is shown in Fig. 13. 

Five different event types can readily be identified in 
this correlation plot. (1) Real three- fold HMS/front- 
array/rear-array coincidence events are denoted "i?" and 
form the peak centered at cTOF = rTOF = ns. (2) 
Three-fold accidental coincidences, denoted "^3" , re- 
quire a random electron in the HMS, a random neutral 
particle in the front array, and a random particle in the 
rear array, and are distributed uniformly over the entire 
plot area. (3) Real two- fold front-array/rear-array co- 
incidences with an accidental electron are denoted ">le" 
and are associated with the "horizontal band" defined 
by rTOF = ns. (4) Real two-fold electron/front-array 
coincidences with an accidental rear array particle are de- 
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FIG. 14: Distributions of the invariant mass W before (cross-hatched) and after (soUd) all cuts except for those on Ap/p, 
Ipmissi, and cTOF at = 1.136 and 1.474 {GeV/cf. The vertical dashed Unes denote the final W < 1.04 GeV/c^ cut. 



noted and are identified with the "vertical band" 

defined by cTOF = ns. (5) Real two- fold clcctron/rear- 
array coincidences with an accidental front-array particle 
are denoted "Ajr" . These events are located along a diag- 
onal band defined (approximately) by cTOF = — rTOF. 
Such events are attributed to the corruption of an oth- 
erwise ii-type event by an accidental front array hit oc- 
curing some time AtA before or after the true interac- 
tion. The values of cTOF and rTOF extracted from 
the data will then be cTOF = cTOFuncorr - Ai^ and 
rTOF = rTOFuncorr + AtA, where the "uncorr" sub- 
script denotes the (true) uncorrupted values. For un- 
corrupted values centered on zero, it then follows that 
cTOF = —rTOF, in accordance with the observed re- 
sult. 



3. Quasielastic event selection 

Real _R-typc coincidence events were selected via tight 
cTOF e [-1, 1] ns and rTOF € [-1,8] ns cuts. As evi- 
dence our cuts selected quasielastic ^H(e,e'n)^H events, 
comparisons of invariant mass spectra, W, obtained be- 
fore and after cuts on Ap/p, |pmiss|, cTOF, and rTOF are 
shown in Fig. 14 for our = i_i36 and 1.474 (GeV/c)^ 
points. After all cuts (except for the additional cut on 
W < 1.04 GeV/c^ itself), these distributions converged 
to fairly narrow peaks centered on the neutron mass. 



4- Extraction of asymmetries from time- of -flight spectra 

One-dimensional projections of cTOF are shown in 
Fig. 15 for our lowest and highest points. Histograms 
of rTOF were accumulated for those events falling within 
the [—1,1] ns peak cTOF window. In addition, his- 
tograms of rTOF were accumulated also for a sampled 
background region of [—8, —2] ns in the cTOF spectrum. 



The signal-to-noise ratios were independent of the state 

of the Charybdis magnet at each of our points. 

Sample rTOF spectra summed over all RU, LU, RD, 
and LD events for cTOF peak events at our lowest and 
highest points are shown in Fig. 16. The asymmetric 
tails on the slow sides are due to scattering from pro- 
tons bound in carbon nuclei and other nuclear reactions, 
and the small satellite peak observed in the = 1.474 
(GeV/c)^ spectrum on the fast side at ~ —2.5 ns is at- 
tributed to quasifree Tr^-production in the scintillators, 
followed by decay and detection of a photon in the rear 
array. Indeed, front-to-rear velocity spectra for these 
events are centered on c. This Tr'^-production peak was 
observed in the = 1.136, 1.169, and 1.474 (GeV/c)^ 
rTOF spectra but was absent in the = 0.447 {GeV/cf 
spectrum, as the energies of those neutrons were below 
threshold. 

The yields for those events falling within the [—1,8] ns 
rTOF window were obtained via peak fitting, with contri- 
butions from the Tr'^-production peak and the flat back- 
ground excluded. These yields were then further cor- 
rected for the contents of the rTOF spectra accumulated 
for the sampled cTOF background region. The desired 
quantities, the physical scattering asymmetries, ^, were 
extracted from the final background-subtracted yields in 
the four decomposed rTOF spectra via the cross-ratio 
technique [71]. In obvious notation, the cross ratio, r, is 
defined to be the ratio of two geometric means. 



NruN 



LD 



NrdNlu 

and is related to the asymmetry ^ via 



1 



r + l VNruNld + \/ NrdNlu ' 



(20) 



(21) 



The merit of the cross-ratio technique is that ^ is insen- 
sitive to [71]: (1) the number of particles incident on the 
polarimeter (i.e., target luminosities) for the two beam 
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FIG. 15: Distributions of cTOF after application of the final set of cuts at = 0.447 and 1.474 (GeV/cf. The dark shaded 
regions indicate the selected peak window, while the cross-hatched regions indicate the sampled background region. 
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FIG. 16: Distributions of rTOF for cTOF peak events at = 0.447 and 1.474 (GeV/c)^. The cross-hatched regions indicate 
the accepted window. The solid curves are the results of our fits to these spectra. 



helicity states; and (2) the relative efBciencies and accep- 
tances of the polarimeter's top and bottom rear arrays. 



G. Asymmetry results 

1. Electron beam polarization normalization 

Unlike recoil polarization measurements in which both 
polarization components, and , can be extracted 
simultaneously from the data (e.g., recoil polarization ex- 
periments with focal-plane polarimeters), om- polarime- 
ter was sensitive to only one of these components (or a 
combination thereof). As such, it was necessary to nor- 
malize our run-by-run scattering asymmetries to some 
common value of the beam polarization. 

Our normalization procedure was as follows. As the 
beam polarization was measured only periodically with 
the M0ller polarimeter, we defined the beam polariza- 



tion for a run to be the result of the most recent prior 
M0ller measurement (if the accelerator parameters were 
unchanged in the interim). All of our run- by-run scat- 
tering asymmetries and their statistical errors were then 
normalized to a common value of 80%. 

We found that the beam polarization was fairly stable, 
with small (few percent) fluctuations observed in suc- 
cessive measurements during periods of continuous beam 
delivery to our experiment. To illustrate, the results of 
23 successive M0ller measurements conducted during our 
= 1.474 (GeV/c)^ x = ±40° running period spanning 
the days of February 20, 2001 through March 5, 2001 are 
shown in Fig. 17. 

2. Corrections for charge-exchange in the lead curtain 

Contamination from the two-step ^H(e, e'p) + Pb(p', n) 

charge-exchange reaction in the lead curtain could either 
dilute the "real" ^H(e, e'n) asymmetry or contribute to a 
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FIG. 17: Results of 23 successive M0ller beam polarization 
measurements conducted during the = 1.474 (GeV/c)'^ 
X = ±40° running period spanning the days of February 20, 
2001 through March 5, 2001. The errors shown are statistical. 



TABLE V: Final (n, n) and (n, p) asymmetry data normalized 
to a beam polarization of 80%. The = 1.474 (GeV/c)^ 
X = 0° asymmetries were corrected for contamination from 
charge-exchange in the lead curtain. 



Central 


Precession 


(n, n) 




[(GcV/c)2] 


Angle X 


C [%] 


e [%] 


0.447 


-40° 


-4.51+0.22 


-2.97 + 0.19 


0.447 


+40° 


6.38 ± 0.28 


4.98 ± 0.29 


1.136 


0° 


1.20 + 0.13 


0.57 + 0.10 


1.136 


-90° 


-5.71 + 0.32 


-3.11 + 0.25 


1.136 


+90° 


5.67 + 0.35 


3.18 + 0.25 


1.169 


-40° 


-2.92 + 0.29 


-1.42 + 0.22 


1.169 


+40° 


4.75 + 0.31 


2.76 ± 0.25 


1.474 


0° 


1.29 + 0.19 


0.64 ± 0.17 


1.474 


-40° 


-2.26 + 0.20 


-0.88 ± 0.18 


1.474 


+40° 


4.03 + 0.24 


2.11 + 0.21 


1.474 


-90° 


-4.64 + 0.47 


-2.92 + 0.50 


1.474 


+90° 


5.07 ± 0.49 


2.14 + 0.43 



false asymmetry if the flux of charge-exchange neutrons 
was unpolarized or polarized, respectively. A significant 
advantage of our neutron flight path setup in which the 
lead curtain was located downstream of the Charybdis 
dipole field was that the majority of the quasielastic pro- 
tons were swept from the front face of the lead curtain. 

Accounting for such nuclear reactions, the measured 
asymmetry, can be parametrized as 

€M = Mfl + /B€B, (22) 

where ]b denotes the contamination level from the two- 
step charge- exchange process, denotes the asymmetry 
for charge-exchange neutrons, /r = 1 — /s denotes the 
fraction of ^H(e, e'n) neutrons, and denotes the asym- 
metry for the ^H(e, e'n) reaction. The asymmetry for the 
background process can further be written as 

= (P|cosxp + P|;sinXp)^^^A„ (23) 

where and P£ denote, respectively, the projections 
of the ^H(e, e'p) recoil proton's polarization on the po- 
larimcter momentum basis S- and L-axis; Xp is the pro- 
ton spin precession angle in the Charybdis field; and 
ly^ss denotes the polarization transfer coefficient for the 
Pb(p, n) reaction. It then follows that if /s, P|, P£, 
Xp, and Ay are all known or measured, can be 

determined. 

To estimate the contamination levels, /b, we took data 
with a hquid hydrogen target. The rates for (n,n) and 
(n,p) events extracted from these data were compared 
with those extracted from our liquid deuterium data and 
corrected for differences in the two targets' densities and 
atomic numbers. We found that the contamination levels 
were negligible (< 0.3%) at all of our Q"^ points when the 
Charybdis field was energized for x = ±40° and ±90° 
precession and also when the field was de-energized at 
= 1.136 (GeV/c)^ for x = 0° precession; therefore, 



we did not apply corrections to any of these asymme- 
tries. Non-negligible event rates were observed when the 
Charybdis field was de-energized for x = 0° precession 
at = 1.474 (GeV/c)^, with observed contamination 
levels of ~ 2.2% and ~ 4.2% for (n,n) and {n,p) events, 
respectively. Corrections were applied to these asymme- 
tries assuming Dgg = for our kinematics of Tp ~ 786 
MeV. D^^ was measured at Tp = 795 MeV and found to 
be consistent with zero (0.014 ± 0.013) [72]. 



3. Summary of asymmetry results 

Our final asymmetry data for (n, n) and {n,p) events at 
each of our Q"^ points and precession angles are tabulated 
in Table V. To illustrate the quality of our asymmetry 
data, a histogram of the (n, n) asymmetries for the = 
1.136 (GeV/c)2 X = 0° data set is shown in Fig. 18; the 
distribution is of an appropriate Gaussian shape. 



H. Extraction of uncorrected values for Gsn/GMn 

We extracted values for Geu/Gmti from our asym- 
metry data assuming elastic scattering from a free neu- 
tron and infinitesimal pointlike HMS and NPOL ac- 
ceptances and neglecting nuclear physics corrections for 
FSI, MEC, and IC. To do so, we fitted the asymmetries 
as a function of the precession angle to the functional 
form ^(x) cx sin(x + S), where the phase-shift parameter 
S = tan~^ (Pj'^^/P^^''^) was defined in terms of form fac- 
tors and kinematics in Eq. (4). To illustrate the quality 
of these fits, our = 1.136/1.169 (GeV/c)^ (n,n) and 
{n,p) asymmetry data are plotted as a function of the 
precession angle in Fig. 19. These data are fitted well 
by sinusoids with excellent agreement seen between the 
independent fits to the (n, n) and (n,p) asymmetry data. 
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TABLE VI: Values of S = tan '^(P/ '''/P^*''') and the uncorrected results for GEn/GMn at each of the points. 



Central 

[{GcV/cf)] 


S [deg 

(n, n) 


(n,p) 


GEn/GMn 

(n, n) 


{n,p) 


GEn/GMn 

Combined" 


0.447 


8.2 ± 1.5 


12.0 ±1.9 


-0.0580 ± 0.0106 


0.0854 ± 0.0138 


-0.0681 ± 0.0084 


1. 136/1. leg*" 


11.7 ±1.2 


11.2 ±1.7 


-0.124 ±0.013 


-0.118 ±0.019 


-0.122 ±0.011 


1.474 


14.0 ± 1.6 


16.9 ±2.9 


-0.166 ± 0.020 


-0.203 ± 0.037 


-0.174 ±0.017 



"Weighted average of Geu/Gmu from (n,n) and (n,p) events. 

'Result obtained via averaging of the nominal (central) electron kinematics for the two points. 



Q^= 1.136 (GeV/cf ; X = 0°; (n,n) 




^ [%] 

FIG. 18: Histogram of the = 1.136 (CeV/c)^ X = 0° 
(n, n) asymmetries. The solid curve is a Gaussian fit, and the 
vertical dashed line denotes the mean value for the asymmetry 
given in Table V. 
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FIG. 19: Sinusoidal fits of the = 1.136/1.169 {GeV/cf 
{n,n) and {n,p) asymmetries as a function of the precession 
angle. 



Wc could not fit the = 0.447 (GcV/c)^ asymmetries 
to a sinusoid as asymmetry data were taken only at two 
precession angles. 

The values for GEn/GMn wc derived from our values 
for 6 using the nominal (central) values for the kinematics 



listed in Table II are summarized in Table VI. 



I. Simulation programs 

We developed two independent simulation programs, 
GENGEN and the Acceptance program, to extract 
acceptance-averaged and nuclear physics corrected values 
for GEn/GMn from our measured experimental asymme- 
tries. The GENGEN simulation program, a pure Monte 
Carlo simulation program, included realistic models for 
the primary ^H(e, e'n)^H reaction in the target, the HMS 
acceptance, neutron spin precession in the Charybdis 
dipole field, spin-dependent neutron scattering in the 
lead curtain, elastic and quasielastic np scattering in the 
front and rear arrays of NPOL, tracking of the incident 
neutron and recoil proton from the front array to the rear 
array, and the detector response of the polarimeter to np 
interactions in the front and rear array. The Acceptaince 
program was not a Monte Carlo simulation program, but 
was, instead, designed to extract the corrections for the 
finite experimental acceptance and nuclear physics effects 
directly from our experimental data. 



1. GENGEN simulation program 

a. Event sampling technique A uniform sampling 
scheme was employed in which events were generated uni- 
formly over the available kinematic phase space, with an 
event weight computed according to a model cross sec- 
tion. The vertex position for the primary ^H(e, e'n)^H in- 
teraction in the extended target was sampled uniformly 
within the raster pattern, and the scattered electron's 
kinematics were sampled uniformly over specified ranges. 
The physical acceptance of the HMS was enforced via in- 
clusion of an HMS transport model taken from the SIMC 
simulation code [73] . In-plane and out-of-plane scattering 
angles for the recoil neutron were sampled uniformly over 
specified ranges, permitting computation of the magni- 
tude of the neutron's momentum according to Eq. (17). 
Complete specification of the electron and neutron kine- 
matics permitted computation of those variables of par- 
ticular interest for the quasielastic ^H(e, e'n)"^H reaction, 
such as 6^™-, pmiss, etc. 
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b. Cross section and recoil polarization We em- 
ployed the Arenhovel formalism [18, 19, 25] for com- 
putation of the ^H(e, e'n)^H differential cross section 
and recoil polarization. These calculations modeled the 
deuteron as a non-relativistic n-p system and employed 
the Bonn i?-space NN potential [74] for the deuteron 
wavefunction and the inclusion of FSI; further, leading- 
order relativistic contributions (RC) to the wavefunc- 
tions and one-body current were added via inclusion of 
the most important kinematic part of the wavefunction 
boost. In the current operator, explicit MEC contribu- 
tions beyond the Siegert operators (essentially from tt- 
and p-exchange) and IC were included. The treatment 
of IC permitted permitted consideration of kinematic re- 
gions away from the quasielastic ridge and excitations up 
to the A region. 

Acceptance-averaging of those calculations performed 
within the Born approximation (hereafter, termed the 
"PWB A model" ) permitted extraction of the corrections 
for the finite experimental acceptance (over the point- 
like results discussed in Section VH), while averaging of 
the full calculations which included FSI, MEC, IC, and 
RC (hereafter, termed the "FSI+MEC+IC+RC model") 
permitted application of corrections for nuclear physics 
effects. In order to implement the Arenhovel formalism 
within GENGEN, lookup tables for the structure functions 
for the ^H(e, e'n)^H reaction were constructed over a suf- 
ficiently dense kinematic grid indexed by {Ef.' , 9e' , ©Ji'™ ), 
and tri-cubic spline interpolation among the grid ele- 
ments was used to compute the cross section and recoil 
polarization for the kinematics of each simulated event 
according to the formalism outlined in Appendix A. 

c. Nucleon form factors All of the structure func- 
tion calculations assumed the dipolc paramctrization for 
Gmti, Gep, and Gup- For the form factor of interest, 
Geu, the structure function calculations were first per- 
formed for various multiplicative factors of the standard 
Galster paramctrization, Gsn = —SiinTGn/i^ + 5.6t), 
where the scale factor S e {0.50,0.75,1.00,1.25,1.50}. 
To investigate the influence of a different depen- 
dence for Geu, structure function calculations were per- 
formed also for multiplicative factors of a modified Gal- 
ster paramctrization, Geu = —SajinTGn/i^ + br), with 
a = 0.894, b = 3.55 (which choice will be explained later), 
and the same set of S factors. 

d. Charybdis field transport The recoil polarization 
was transported point-by-point through a grid of the 
Charybdis field, with the time derivative of the spin vec- 
tor computed at each grid point according to standard 
relativistic electrodynamics. The precession angle was 
computed from information on the initial and final spa- 
tial orientations of the spin vector. 

e. Lead curtain interactions Neutron interactions in 
the lead curtain were simulated with a spin-dependent 
multiple scattering algorithm that employed quasifree 
scattering from a lead nucleus modeled as a Fermi gas, 
with the Fermi momentum for ^"^Pb taken to be 265 
MeV/c [2]. The probability for an interaction of the neu- 



tron with a lead nucleus was determined via interpola- 
tion (or extrapolation) of existing data on total n + Pb 
cross sections [75] . A polar scattering angle was sampled 
from cumTilative probability distributions for the polar 
scattering angle as a function of neutron energy, and an 
azimuthal scattering angle was chosen via an acceptance- 
rejection algorithm for the spatial scattering asymmetry 
resulting from non-zero analyzing power. Pauli blocking 
was enforced. For those neutrons suffering an interac- 
tion, the scattered neutron's and recoil nucleon's polar- 
ization components were constructed via computation of 
the depolarization and polarization-transfer tensors for 
NN scattering using helicity amplitude routines obtained 
from SAID [70]. 

/. Polarimeter interactions Finally, following (suc- 
cessful) transport of the neutron through the steel col- 
limator into the front array, interactions in NPOL were 
simulated. A scattering vertex was chosen randomly as- 
suming a fixed value for the mean free path of neutrons 
in the plastic scintillator, and both the clastic (scatter- 
ing from free protons) and quasielastic (scattering from 
protons bound in carbon nuclei) channels were simulated. 
The scattering angles in the polarimeter were determined 
using the same algorithms employed for NN scattering in 
the lead curtain. We employed a rather simple model for 
the propagation of the recoil proton, with the energy de- 
position and range (assuming a straight-line trajectory) 
computed according to the Cecil, Anderson, and Madey 
[76] range-energy formulas for protons in the hydrocar- 
bon scintillator. 



2. GENGEN performance 

A rigorous and reliable extraction of the corrections for 
the finite experimental acceptance and nuclear physics 
effects from simulated data is feasible if the simulated 
acceptance reasonably matches the experimental accep- 
tance; therefore, we now document the performance of 
GENGEN by comparing: (1) simulated distributions of im- 
portant kinematic quantities with those derived from ex- 
perimental data; and (2) the behavior of the acceptance- 
averaged simulated polarizations and the experimental 
asymmetries as a function of the cut on some kinematic 
variable (here, taken to be the invariant mass W). 

a. Kinematic distributions Sample comparisons of 
experimental and simulated kinematic distributions of 
two important kinematic variables, W and [pmissj) are 
shown in Figs. 20 and 21. Reasonable agreement is seen 
between the GENGEN distributions and those extracted 
from experimental data. Although not shown here, rea- 
sonable agreement was also obtained between simulated 
and experimental distributions of variables related to np 
scattering in NPOL (e.g., scattering angles, velocity spec- 
tra, etc.). 

b. Experimental asymmetries and simulated polariza- 
tions A sample comparison of the behavior of the exper- 
imental asymmetries and acceptance-averaged simulated 
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FIG. 20: Comparison of GENGEN simulated (unfilled histograms with thick solid line borders) and experimental (cross-hatched 
filled histograms) distributions of W for the four central points. Identical cuts were applied to both the simulated and 
experimental data. The simulated results shown here employed the FSI+MEC+IC+RC model and the Galster parametrization 
for Geu- 



polarizations following transport through the Charybdis 
dipole field is shown in Fig. 22. There, we plot the ratio of 
the experimental asymmetries to the simulated polariza- 
tions as a function of the upper cut on W for {n, p) events 
and X = -40° precession at our Q"^ = 1.169 (GeV/c)^ 
point. Within statistical errors, the experimental asym- 
metries and simulated polarizations arc seen to scale sim- 
ilarly with the cut on W . Similar results were observed 
for our other Q'^ points and precession angles. 



It should be noted that in this figure the simu- 
lated acceptance-averaged polarizations were computed 
assuming some certain parametrization for Geu (hero, 
the Galster parametrization); therefore, the ratios of the 
asymmetries to the simulated polarizations shown in this 
figure are not equivalent to the polarimeter's analyzing 
power. 



3. Acceptance program 

The Acceptance program was developed as an alter- 
native to the GENGEN Monte Carlo simulation program. 
This program used the kinematics of the reconstructed 
quasielastic events from the actual experimental data to 
compute, on an event- by- event basis, the recoil polariza- 
tion presented to the polarimetcr for each event employed 
in our final data analysis (i.e., for those events satisfying 
all final analysis cuts). The Acceptance program used 
the same ^H(e, e'n)^H interpolation and Charybdis spin 
transport algorithms developed for GENGEN. Although the 
Acceptcince program was, technically, not a true Monte 
Carlo simulation, a significant advantage of this method 
was that it did not require a model for the experimen- 
tal acceptance; however, the disadvantage of this method 
was that the reconstruction of the event- by- event kine- 
matics is, of course, subject to measurement uncertain- 
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FIG. 21: Comparison of GENGEN simulated (unfilled histograms with thick solid line borders) and experimental (cross-hatched 
filled histograms) distributions of |pmiss| for the four central points. Identical cuts were applied to both the simulated and 
experimental data. The simulated results shown here employed the FSI+MEC+IC+RC model and the Galster parametrization 
for Geti- 



ties, leading to uncertainties in the computation of the 
recoil polarization. 



VI. FINAL RESULTS FOR Gs-alGM-a AND Geu 

A. Distributions of ©^■"• 

Distributions of 6$^^' for those events surviving all 
analysis cuts at our lowest point are shown in Fig. 
23. The majority of the accepted events are seen to fall 
within ~ 10-15° of perfect quasifree emission. The dis- 
tributions of 6^^' at our other points are similar, 
but arc restricted to somewhat smaller ranges, 170° < 

e^^-- < 180°. 

Even for perfect quasifree emission, O^™' = 180°, 
the PWBA and FSI+MEC+IC+RC calculations of the 
jp^^ polarization ratio differ by 4.2% for the central 



kinematics of our lowest = 0.447 (GcV/c)^ point and 
1.6% at our highest = 1.474 (GeV/c)^ point. As the 
differences between the PWBA and FSI+MEC+IC+RC 
calculations increase away from 0^'™' = 180°, these num- 
bers provide essentially lower bounds for the expected 
magnitude of corrections for nuclear physics effects. 



B. Extraction of acceptance-averaged and nuclear 
physics corrected values for GEn/GMn 

1. Overview of acceptance- averaging analysis procedure 

The recoil polarization component we were interested 
in was the projection of the polarization vector on the 
polarimeter momentum basis S-axis following transport 
through the Charybdis field and the lead curtain. We de- 
note this polarization component as Pg, where the prime 
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FIG. 22: Ratio of the asymmetries extracted from the ex- 
perimental data to the GENGEN simulated polarizations as a 
function of the cut on W for (n,p) events and x = —40° 
precession at our = 1.169 (GeV/c)^ point. The shaded 
band indicates the statistical error on the ratio for the nomi- 
nal cut on of < 1.04 GeV/c^. The simulated results shown 
here employed the FSH-MEC+IC+RC model and the Galster 
parametrization for Gb„. 



q2 = 0.447 {GeV/cf 



3000 



2000 



1000 




166 168 170 172 174 176 178 180 

jC.m. 
np 



©"'"■ [deg] 



FIG. 23: Distributions of OJip" after application of the final 
set of analysis cuts at = 0.447 (GeV/c)^. 



and X = ±40°) and then computed a value for each 
precession angle combination and each event type [i.e., 
(n, n) or (n, p) events] according to 



^cxp ) 



(A?7sim) + (A?7cxp) 



(24) 



where r/sim = (P^(0°))/(P^(±90°)) for the 
X = 0°,±90° precession angle combination and 
(P^(-40°))/(P^(+40°)) for the X = ±40° precession an- 
gle combination. The expressions for Tjexp are identical, 
with the acceptance-averaged polarizations replaced by 
the experimental asymmetries. Arysim and A?7cxp denote 
the statistical errors. The resulting x^ values were 
fitted as a function of the scale factor S to a parabolic 
function, with the optimal value of S defined by the zero 
of the parabolic fitting function. 

(3) In our "global analysis method" , we compared the 
experimental asymmetries with the simiilatcd polariza- 
tions via minimization of a global value computed 
according to 



E 



A 



{n,n), (n,p) 



(AO^ + (A(P^))^ 



(25) 



Here, the sum runs over all 10 asymmetries, ^, and sim- 
ulated polarizations, (Pg), for each point [i.e., 5 dif- 
ferent precession angles, and (n, n) and (n,p) events], 
and A^'^^ and A^'^^ denote the polarimetcr's analyzing 
power for (n, n) and {n,p) events. A^ and A(P^) de- 
note the statistical errors. In this analysis, the analyzing 
powers and scale factor S were treated as free parame- 
ters, with the optimal values extracted from the minimal 
value. 

We note that the simulation statistical errors were gen- 
erally an order of magnitude smaller than the experimen- 
tal statistical errors. 



denotes transport through the dipole field and lead cur- 
tain. Acceptance-averaged and nuclear physics corrected 
values for GEn/Gnin were extracted from our experimen- 
tal asymmetries and simulations at each point via the 
following procedure: 

(1) Acceptance-averaged polarizations (Pg) computed 
according to the PWBA and FSI+MEC+IC+RC models 
were extracted from simulated data for each precession 
angle at each point and for each scale factor S of the 
Galster parametrization (see Section VII). 

(2) In our "pairwise analysis method", for each S fac- 
tor, we compared the ratio of the experimental asymme- 
tries to the ratio of the simulated polarizations for the dif- 
ferent precession angle combinations (i.e., x = 0°,±90° 



2. Acceptance- averaged values of 

The acceptance-averaged values of Q^, denoted (Q^), 
were determined to be (Q^) = 0.447, 1.126, 1.158, and 
1.450 (GcV/c)2 for the central = 0.447, 1.136, 1.169, 
and 1.474 (GeV/c)^ points, respectively. The distribu- 
tion of Q2 values for the (Q^) = 0.447 {GeV/cf point 
was sharply peaked around the central value of 0.447 
(GeV/c)^, while the distributions of values for the 
(Q2) = 1.126/1.158 and 1.450 (GeV/c)^ points were in- 
tegrated from ^ 1.0 to ^ 1.3 (GeV/c)^ and from ^ 1.2 
to ~ 1.7 (GeV/c)^, respectively. 

Henceforth, we will use (Q^) = 1.132 (GeV/c)^ to 
denote the sample-size weighted average of the (Q^) = 
1.126 and 1.158 (GeV/c)^ data sets. 
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3. Acceptance-averaging analysis iterations 

We performed two iterations of the above-described 
analysis procedure with both the Acceptance and 
GENGEN simulation programs. 

In the first iteration, the simulations were conducted 
with the PWBA and FSI+MEC+IC+RC calculations 
that assumed different multiplicative factors of the stan- 
dard Galster parametrization for the dependence of 
Geu- The optimal values for the scale factors S were 
then used to compute the optimal values for GeuIGmu 
according to GEnlGun = -S'optimai X (t)/(1 + 5.6(t)), 
where (r) = ((5^)/4rn^. Values for GEn were then ex- 
tracted from our optimal values for Geu/Gmu using the 
best-fit values for Gmk taken from the parametrization 
of Kelly [49] . Then we fitted our first-iteration results for 
Geu together with the then-available world data on Geu 
(as of early 2003) to the modified Galster parametrization 
described previously in Section VII; the best-fit param- 
eters we found at that time were a = 0.894 ± 0.023 and 
b = 3.55 ± 0.37. This fit included the then-available data 
on Geu extracted from measurements using polarization 
degrees of freedom [33, 38, 40, 41, 43, 48] and an analysis 
of the deuteron quadrupole form factor [30] , and also data 
on the slope of Geu as measured via low-energy neutron 
scattering from electrons in heavy atoms [77]. Since this 
analysis, new data on Geu have been published [50, 51], 
and a new modified Galster parametrization has been 
published [78]. 

In our second analysis iteration, a second set of the 
PWBA and FSI-|-MEC-hIC-l-RC calculations were per- 
formed that assumed this modified Galster parametriza- 
tion for the dependence of G_e„. The Accepteince 
and GENGEN simulations were both repeated using these 
new calculations, and the procedure for the extraction of 
the optimal Geu/Gmti values was identical to that of the 
first iteration. 

The differences between the first and second analysis 
iterations were negligible. This result is not surprising, 
because: (1) both parametrizations have small second 
derivatives in the vicinity of our points; and (2) the 
acceptance was fairly symmetric about the acceptance- 
averaged values of Q"^. 

4- Acceptance-averaging analysis results 

The pairwise analysis method was employed for the 
extraction of our Geu/Gmti values at (Q^) = 0.447 
(GeV/c)^ (only two precession angles), while the global 
analysis method was employed for the analysis of our 
(0^) = 1.132 and 1.450 (GeV/c)^ data sets. The final 
acceptance-averaged and nuclear physics corrected values 
for Geu/Gmii we obtained with the Acceptance program 
and GENGEN agreed to better than 1% at (Q^) = 0.447 
and 1.132 (GeV/c)^ and 2% at (Q^) = 1.450 (GeV/c)^, 
well within the statistical errors; therefore, the values 
for GEn/GMn we report later in Table IX are the aver- 



TABLE VII: Analyzing powers for (n, n) and (n, p) events at 
each of our points. The errors are statistical. 



Event 




{Q') [(GeV/c)^] 




Type 


0.447 


1.132 


1.450 


(n, n) 


0.141 ±0.004 


0.137 ±0.010 


0.144 ±0.013 


{n,p) 


0.103 ± 0.005 


0.075 ± 0.007 


0.071 ±0.011 



age of the central values obtained with our two simula- 
tion programs. The analyzing powers we extracted from 
our acceptance-averaging analysis procedures are sum- 
marized in Table VII. 



C. Systematic uncertainties 

An itemized summary of estimates for the magnitudes 
of our relative systematic uncertainties in Geu/Gmti ap- 
pears in Table VIII. Our final values for the total relative 
systematic uncertainties, 2-3%, are much smaller than 
our relative statistical uncertainties. Brief discussions of 
each itemized systematic uncertainty (and others deemed 
negligibly small) appear below. 



1. Beam polarization 

The beam polarization cancels in the form factor ratio 
only if it does not vary during sequential measurements of 
the scattering asymmetries. Consequently, fluctuations 
in the beam polarization measurements introduce a sys- 
tematic uncertainty. We estimated the temporal uncer- 
tainty in the beam polarization via the following proce- 
dure. First, polarization measurements conducted under 
similar conditions at the polarized source were grouped 
into clusters. Second, the mean value of the polarization 
for each cluster was computed and then recentered about 
the nominal 80% polarization. Next, the statistical error 
for the entire data set (i.e., all identified clusters) was 
computed, and the overall uncertainty was then increased 
by the square root of (to account for the observed fluc- 
tuations). Finally, our total estimated uncertainty in the 
polarization was propagated through the expression for 
the form factor ratio, Eq. (4). 



2. Charge-exchange in the lead curtain 

Estimates of the contamination levels from the two- 
step ^H(e, e'p) + Pb(p, n) charge-exchange reaction were 
given previously in Section VG2. To estimate the sys- 
tematic uncertainty in GEnlGun due to contamination 
from this background process, we computed values for 
the recoil proton's polarization using values for G^p and 
Gmp taken from the parametrization of [79]. These po- 
larization components were then transported through 
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TABLE VIII: Compilation of our estimated relative system- 
atic uncertainties in Geji/Gmti [%]■ The total systematic 
error that is quoted for each point and precession angle 
combination is the quadrature sum of the itemized systematic 
uncertainties. 



m [{GeY/cf] 



Source 


0.447° 


1.132° 1.132' 


1.450" 


1.450'' 


Beam polarization 


1.6 


0.7 


0.4 


1.2 


0.3 


Charge-exchange 


<0.1 


<0.1 


0.1 


<0.1 


0.2 


Depolarization 


<0.1 


0.1 


<0.1 


<0.1 


0.6 


Positioning/traceback 


0.2 


0.3 


0.3 


0.4 


0.4 


Precession angle 


1.1 


0.3 


0.1 


0.5 


0.1 


Radiative corrections 


0.7 


0.1 


0.1 


0.1 


0.1 


Timing calibration 


2.0 


2.0 


2.0 


2.0 


2.0 


Total 


2.9 


2.2 


2.1 


2.4 


2.2 



"'X = ±40° precession. 
= 0°, ±90° precession. 



the Charybdis dipole field using estimates for the pro- 
ton spin precession angles. As there are very few data 
on the lead polarization transfer coefficient, -D55, we 
calculated the correction to the asymmetries (using in- 
formation on the analyzing powers extracted from our 
acceptance-averaging analysis and the values for, and the 
uncertainties in, the charge-exchange contamination lev- 
els) for various (reasonable) choices of -Dg^. Spreads in 
the resulting values of Gsn/GMn were then defined to 
be the systematic uncertainties. 

3. Neutron depolarization in the lead curtain 

The total n -|- Pb cross section is fairly flat at ~ 3 
barns over the range of neutron kinetic energies in our 
experiment (slow rise with energy) [75]. For our 10.16- 
cm thick lead curtain, our GENGEN simulations indicated 
a 30.8%, 42.5%, 43.0%, and 46.7% interaction probabil- 
ity for the neutron energies at our (Q^) = 0.447, 1.126, 
1.158, and 1.450 (GeV/c)^ points, respectively. We found 
that the contamination levels within our [—1, 1] ns cTOF 
window from neutrons suffering one or more interactions 
in the lead curtain were 0.04%, 3.8%, 4.2%, and 9.3% at 
(0^) = 0.447, 1.126, 1.158, and 1.450 (GeV/c)^, respec- 
tively. The fact that our simulations predicted a much 
more rapid increase in the contamination levels with en- 
ergy as compared to the interaction probabilities is be- 
cause the angular distributions for nn and np scattering 
peak at large (small) scattering angles for the neutron 
kinetic energies at our lowest (highest) point (as com- 
puted by SAID [70]). Further, our simulations suggested 
that interactions in the lead curtain may have been partly 
responsible for the small tail observed on the slow side of 
our experimental cTOF distributions at our highest 
point (see Fig. 15). 

The quantity of interest was the spectrum of the po- 
larization presented to the polarimetcr front array for 
neutrons which did and did not interact with the lead 



(Q^) = 1.450 (GeV/c)^ 




P' (x = ^0°) 

o 

FIG. 24: Sample GENGEN simulated (Pg) spectrum for x = 
-40° precession at (Q^) = 1.450 {GeV/cf. The unfilled his- 
togram is summed over all simulated events, while the cross- 
hatched histogram is summed over those events suffering one 
or more interactions in the lead curtain. The units of the 
ordinate are arbitrary. 



curtain. A sample result comparing polarization spectra 
for these two types of events for x = —40° precession at 
(Q^) = 1.450 (GeV/c)2 is shown in Fig. 24. Our sim- 
ulations indicated that the distribution of polarizations 
for neutrons suffering an interaction in the lead curtain 
is a broad continuum, yielding a depolarization of the 
neutron flux presented to the polarimeter. Similar re- 
sults were observed at our other points. We found, 
though, that the effects of depolarization in the lead cur- 
tain tend to cancel in the polarization ratio, leading to 
small systematic uncertainties in the Gsn/GMn ratio. 
The magnitudes of the residual non-cancellations were 
taken to be the uncertainties listed in Table VIII. 



4- Positioning and traceback 

Two contributions to an uncertainty in the electron 
scattering angle were considered: positioning (offset in 
the scattering angle from the nominal value) and trace- 
back (reconstruction from the focal plane to the target). 
For the purposes of this analysis, we assumed the un- 
certainties in the electron scattering angle, AOf.', were 
Adf,/ = 1.2 mrad and 1.3 mrad for the positioning and 
traceback uncertainties, respectively; these values were 
derived from a systematic analysis of kinematic data 
taken during this experiment. The systematic uncertain- 
ties in Geti/Gmu were obtained via propagation of these 
values for A9e' through Eq. (4) for the form factor ratio. 
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TABLE IX: Summary of our final results for GEn/GMn and Geu- The first (second) set of errors is statistical (systematic). 
The results reported here are the weighted average of (n, n) and (n,p) events in the polarimcter. 

(Q^) [(GeV/c)^] 

Analysis Quantity 0.447 1.132 1.450 

n(e, e'n) GEn/GMn -0.0681 ± 0.0084 ± 0.0020 -0.122 ± 0.011 ± 0.003 -0.174 ± 0.017 ± 0.004 

^H(e, e'nfn PWBA GEn/Gun -0.0713 ± 0.0086 ± 0.0021 -0.126 ± 0.010 ± 0.003 -0.183 ± 0.018 ± 0.004 

^H(e,e'n)^H FSI+MEC+IC+RC Geu/Gmu -0.0755 ± 0.0089 ± 0.0022 -0.131 ± 0.011 ± 0.003 -0.189 ± 0.018 ± 0.004 

Values from [78] GMu/jJ-nGp 1.003 ± 0.005 1.067 ± 0.012 1.064 ± 0.016 

n(e, e'n) G^ 0.0492 ± 0.0061 ± 0.0015 0.0370 ± 0.0032 ± 0.0009 0.0383 ± 0.0038 ± 0.0011 

2H(e,e'n)iH PWBA Geu 0.0515 ± 0.0062 ± 0.0015 0.0381 ± 0.0032 ± 0.0009 0.0403 ± 0.0039 ± 0.0011 

2H(e, e'n)iH FSI+MEC+IC+RC Geu 0.0545 ± 0.0064 ± 0.0016 0.0396 ± 0.0032 ± 0.0010 0.0415 ± 0.0039 ± 0.0011 



5. Precession angle 

Uncertainties in the neutron spin precession angle were 
estimated [64] via a calculational scheme that employed 
the reconstructed kinematics from the experimental data 
as the source of the neutron momentum vectors incident 
on the Charybdis dipolc field. Spin vectors were trans- 
ported through the field using the same magnetic spin 
transport algorithms developed for our two simulation 
programs. This technique provided a measure of the sen- 
sitivity of the precession angle to details of the field map. 
The uncertainties in the mean values of the precession 
angles derived from these studies (at the level of ±0.2°) 
were combined in quadrature with two other sources of 
uncertainty. First, as discussed in Section IV C, we ob- 
served small differences between the measured field in- 
tegrals for opposite magnet polarities and also between 
the field integrals derived from our measured maps and 
the calculated TOSCA maps. These uncertainties were es- 
timated to bo on the level of ±0.3°. Second, as also dis- 
cussed in Section IV C, the field was mapped only along 
the central axis; therefore, we assigned further uncertain- 
ties (at the level of ±0.2°) for incomplete knowledge of 
the field beyond the central axis. Our best estimates of 
the total uncertainties in the precession angle were then 
propagated through the form factor ratio, Eq. (4). 



6. Radiative corrections 

Radiative corrections were calculated specifically for 
the kinematics of this experiment by Afanasev et al. [80] . 
The primary effect of radiative corrections on the recoil 
polarization components P/''^ and Pf^^ was found to be 
depolarization of the electron such that both components 
of the recoil polarization should be increased by ~ 1.9%, 
~ 3.7%, and ~ 4.4% at (Q^) = 0.447, 1.132, and 1.450 
(GeV/c)^, respectively; however, these corrections nearly 
cancel in the form factor ratio such that the net effect is 
small at (Q^) = 0.447 (GeV/c)^ and negligible at the two 
higher points. The residual non-cancellations of the 
corrections in the form factor ratio were taken to be the 
systematic uncertainties we quote in Table VIII. 



7. Timing calibration of the polarimcter 

The timing calibrations we deemed suitable for cer- 
tain running conditions (e.g., periods in between changes 
to the high-voltages for the PMTs) were obtained using 
a subset of the; data for that particular running period. 
To assess the dependence of our results for the scattering 
asymmetries on the choice of the subset of data employed 
for the timing calibration, various calibrations were gen- 
erated from difi'erent subsets of the available data. Excel- 
lent agreement was always found between the results for 
the scattering asymmetries obtained from analyses using 
these different calibrations; however, we did find a ~ 2% 
sensitivity of our results to the choice of the subset of 
data employed for the timing calibration. 

8. Other uncertainties 

We deemed two other possible sources of systematic 
uncertainties to be negligible. First, we demonstrated 
quantitatively that our scattering asymmetries were in- 
sensitive (within statistical errors) to a possible geomet- 
ric asymmetry in the polarimeter (i.e., a spin-averaged 
"top-bottom" asymmetry) by varying our software en- 
ergy thresholds on the top (bottom) rear array while 
maintaining a constant threshold on the bottom (top) 
rear array. Second, analysis of our data taken with the 
"dummy targets" (see Section III E) showed that the level 
of contamination within our [—1,1] ns cTOF window 
from scattering in the target cell windows was negligi- 
ble (< 0.05%). 

D. Summary of final GEn/GMn and Gsn results 

Our final results for GEn/GMn and Geu extracted 
from three different analyses are tabulated in Table 
IX and compared in Fig. 25. The three analyses are 
for: (1) clastic n(e, e'n) scattering and infinitesimal 
HMS and NPOL point acceptances; (2) quasielastic 
^H(e, e'n)^H scattering and acceptance-averaging of the 
PWBA model; and (3) quasielastic ^II(e, e'n)'^H scatter- 
ing and acceptance-averaging of the FSI±MEC±IC±RC 
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* n(e,en) 

• ^H(e,e'n): PWBA 
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FIG. 25: Comparison of our results for Geu at (Q^) = 0.447, 
1.132, and 1.450 (GeV/c)^ extracted from the various analyses 
summarized in Table IX. The data points shown for the three 
analyses at each (Q^) point have been slightly displaced about 
the actual (Q^) value for clarity. The solid curve is the Galster 
parametrization [36]. 

model. 

We note that in our first publication [11] we used val- 
ues for Gun taken from [49]; here we use slightly differ- 
ent values for Gun taken from [78] . The total systematic 
errors we quote for Geu are the quadrature sum of the 
experimental sytematic errors and the relative uncertain- 
ties in Gmu- 

E. Comparison of nucleon form factor data with 
selected theoretical model calculations 

The availability of precise new data on nucleon form 
factors has stimulated much more theoretical work in 
the last few years than we can review here; our se- 
lection of models is not intended to be complete. Al- 
though the present experiment is limited to G^n, we 
believe that comparison with models must consider all 
four Sachs form factors. In this section, we compare rep- 
resentative models with selected data. The data from 
this experiment are shown as filled circles in Figs. 26-29 
while data from other experiments are shown as open cir- 
cles. We selected [81, 82] for Gsp, [83-86] for Gep/Gmp, 
[47, 79] for Gmp, [38, 40, 41, 43, 48, 50, 51] for Gb„, and 
[8, 44, 46, 87, 88] for Gmu- 

1. Models based upon vector meson dominance 

Models based upon vector meson dominance (VMD) 
postulate that the virtual photon couples either directly 

to an intrinsic nucleon core or through the mediation of 
a neutral vector meson, initially limited to the lowest w, 
p, and (j) mesons. The earliest versions assumed that the 
core is a structureless Dirac particle. lachcllo et al. [89] 
assigned the core a form factor and employed a model of 



the p width. Gari and Kriimpelmann [90, 91] then per- 
mitted different Dirac and Pauli intrinsic form factors 
and introduced modifications that ensured consistency 
with pQCD scaling at large and logarithmic running 
of the strong coupling constant. Bijker and lachello [92] 
adopted the Gari and Kriimpelmann (GK) pQCD pre- 
scriptions and refit their model to modern data, still us- 
ing a common intrinsic form factor. This fit, using a 
total of 6 free parameters, is compared with the data in 
Fig. 26. Finally, Lomon [93, 94] produced a more flexible 
set of fits using a model described as "GK extended"; 
the GKcx(02S) version is also shown in Fig. 26. The 
Lomon model uses two intrinsic form factors, the GK 
prescription for the pQCD Hmit, and includes p'(1450) 
and a;'(1419) couplings in addition to the customary p, 
w, and (j) couplings. The p width is included but the 
p' and ijj' structures are not. The fit achieved by this 
extended model, with 13 free parameters, is clearly su- 
perior, especially at large Q^. The Bijker and lachello 
model describes the qualitative behavior of Gsp, but its 
transition between Gep/Gd ~ 1 at low and the nearly 
linear decrease for 1 < < 6 (GeV/c)^ is too grad- 
ual. Nor docs it reproduce the slope in Gupl P-pGu for 
Q"^ > 10 (GeV/c)^. Both of these features are fit well 
by the Lomon model. Unfortunately, the neutron data 
do not discriminate between these models very strongly. 
The Bijker and lachello model provides a slightly bet- 
ter fit to the present Geu data, but the Lomon fit was 
performed before these data became available; it is likely 
that only a slight parameter adjustment would be needed 
to achieve a comparable fit without sacrificing the fits to 
other form factors. It will be interesting to see whether 
the rather large values for Geti/Gd for Q'^ > 2 (GeV/c)^ 
predicted by the Bijker and lachello model are confirmed 
by upcoming experiments [95, 96] that will probe Geu to 
= 4.3 {GeV/cf. Note, however, that the Bijker and 
lachello fit is systematically above the Gmu data for the 
same kinematics, > 2 (GeV/c)^. 



2. Models emphasizing the pion cloud 

The role of the pion in mediation of the long-range 
nucleon-nucleon interaction clearly demonstrates its im- 
portance in understanding form factors for low Q^. Typi- 
cal pion cloud models describe nucleon form factors using 
diagrams in which the virtual photon couples to either a 
bare nucleon core or to the nucleon or the pion loop in 
a single-pion loop. Some models also permit excitation 
of the intermediate state and include additional contact 
terms. A relatively simple example is the Adelaide ver- 
sion [97] of the cloudy bag model (CBM) in which the 
core is based upon the bag model, intermediate excita- 
tion is neglected, and relativistic corrections are made 
using a simple ansatz for Lorentz contraction [98]. Pre- 
dictions from Lu et al. [97] using a bag radius of 0.8 
fm are compared with the data in Fig. 27. Although 
density-dependent extensions of this model, described as 
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FIG. 26: Comparison of representative VMD models with nucleon form factor data {Gsp from [81, 82]; Gep/Gm-p from [83-86]; 
Gmp from [47, 79]; Ge^ from [38, 40, 41, 43, 48, 50, 51]; Gmu from [8, 44, 46, 87, 88]). Dashed curve: Bijker and lachello [92]. 
Solid curve: version GKex(02S) of Lomon [94]. 
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FIG. 27: Comparison of representative pion cloud models with nucleon form factor data [Gsp from [81, 82]; Gep/Gmp from 
[83-86]; Gmp from [47, 79]; Ge„ from [38, 40, 41, 43, 48, 50, 51]; Gmu from [8, 44, 46, 87, 88]). Short-dashed curve: QMC 
model [97]. Solid curve: LFCBM [103]. Long-dashed curve: Friedrich and Walcher parametrization [52]. 
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the quark-meson coupling (QMC) model, have been used 
to study the sensitivity of recoil polarization in nucleon 
electromagnetic knockout to medium modifications of the 
nucleon form factors [99-102], its description of free form 
factors is rather poor and one must hope that the density 
dependence of Ge/Gm ratios is more accurate. 

Alternatively, the light front cloudy bag model 
(LFCBM) of Miller [103] maintains Poincare invariance 
by formulating wave func;tions using the light-front ap- 
proach. This version should then be applicable to higher 

. There are only 4 adjustable parameters and the re- 
sults for Set 1 are compared with data in Fig. 27. A 
previous version of this model [104] provided one of the 
earliest predictions of the sharp slope in Gsp/Gup for 

> 1 (GcV/c)^, but the agreement with recent recoil- 
polarization data is only qualitative. The LFCBM cal- 
culation for Gmp/Gd also decreases too rapidly at large 
Q^. Calculations using this model agree relatively well 
with the Geu data for > 1 (GeV/c)^ but are too 
small at lower . Interestingly, this model predicts 
much stronger values for Geu/Gd at large than the 
Lomon parametrization. However, the LFCBM calcula- 
tions for three of the four form factors show complicated 
and rather implausible shapes for < 1 (GeV/c)^ that 
disagree strongly with data. 

Chiral effective field theory [105] provides a more sys- 
tematic procedure that includes intermediate excitation 
and can be extended to two pion loops [106]. Alterna- 
tively, two-loop contributions can be evaluated in disper- 
sion theory [107]. Recently it has become possible also to 
include both pion loops and vector meson diagrams in a 
consistent manner [108]; however, we do not show curves 
here because this approach remains limited to < 0.4 
(GeV/c)2. 

Friedrich and Walcher [52] performed a phenomenolog- 
ical analysis of the nucleon electromagnetic form factors 
using a parametrization motivated by pion cloud models. 
The core form factor is represented by two dipole form 
factors with different ranges while the pion cloud contri- 
bution, represented as a "bump" at low Q^, is described 
by two Gaussians. These fits, with 5 free parameters for 
Geu and 6 for each of the other form factors, are also 
compared with data in Fig. 27. The quality of these fits 
is generally satisfactory, but it is not clear that the pos- 
tulated oscillation in Geu/Gd is warranted by the avail- 
able data; considerably better experimental precision at 
~ 0.3 (GeV/c)^ would be needed to justify such a 
structure. 

A closer look at the Geu data is given in Fig. 28. The 
original Friedrich and Walcher fit (short-dashed curve) 
used a very preliminary version of the data from the 
present experiment and falls systematically below the fi- 
nal data for this and other more recent experiments for 

> 0.5 (GeV/c)^. A re-analysis using final data for 
this experiment plus new data [43, 50, 51] was made by 
Glazier et al. [51] and is shown as the long-dashed curve 
featuring a bump for ^ 0.3 (GeV/c)^ superimposed 
upon a much flatter core form factor. With 5 parameters 




0.5 



1.0 1.5 
[(GeV/c)^] 



2.0 



FIG. 28: Closer look at comparison of representative pion 
cloud models with data on Gsn (data from [38, 40, 41, 43, 
48, 50, 51]). Solid curve: a fit based upon the pion cloud 
model of Kaskulov and Grabmayr [109] . Short-dashed curve: 
parametrization of Friedrich and Walcher [52] . Long-dashed 
curve: re-analysis by Glazier et al. [51] using the Friedrich and 
Walcher model. The dash-dotted curve is the original Galster 
parametrization [36]. 



it is obviously possible to fit the data very well, perhaps 

too well the simple two-parameter fit of Kelly [78] 
based upon the Galster parametrization already provides 
xt = 0.8 without distinguishing between soft and hard 
structures. The data presently available do not require 
this complication. Data at higher should test whether 
such a hard core is needed but significantly more precise 
data for low would be needed to establish the soft 
pion cloud contribution to Geu- 

Finally, Kaskulov and Grabmayr [109] used a chiral 
quark model (xQM) to derive a relationship 



Geti = S{1 - F.^)Gc, 



(26) 



between Geu, the pion form factor F-^, and the core 
form factor Gc for the 3-quark component of the nu- 
cleon. The coefficient 5 is a weighted average over spec- 
troscopic factors for N and A intermediate states in 
the one-pion loop contribution to the self-energy but is 
treated as an adjustable parameter. If one stipulates 
a monopole for F^r = (1 -l- Q^/A^)"^ and a dipole for 
Gc = (1 + Q'^/^c)~^' the neutron electric form factor 



Geu — S 



br 



l + br 



Gc, 



(27) 



with b — 4m^ / A^ reduces to a Galster-like form with up 
to 3 free parameters (5, A^, Ac); however, S is largely 
determined by the neutron radius 



dG 



En 



ag2 



iSb 
2m%' 



(28) 



If we further assume that A^^ within a loop is the same as 
that for pion electroproduction, only Ac remains to be 
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FIG. 29: Comparison of representative quark models with nucleon form factor data {Gep from [81, 82]; Gep/Gmp from [83 86]; 
Gmp from [47, 79]; Gsn from [38, 40, 41, 43, 48, 50, 51]; Gmu from [8, 44, 46, 87, 88]). Solid curve: PFSA using pointlike 
constituents [110]. Long-dashed curve: light-front using OGE interaction and constituent-quark form factors [111]. Dash-dotted 
curve: hCQM with constituent-quark form factors [113]. 



fit to data for Geu- Thus, using fixed parameters S = 
0.26 and b = 6.65 suggested by Kaskulov and Grabmayr, 
wc fit = 1.00 ± 0.03 (GcV/c)2 to the current Geu 
data. The value given in [109] for Ac is slightly smaller 
because they used the same preliminary data as [52] that 
are smaller than the final results. Our fit is shown in 
Fig. 28 and is practically indistinguishable from the two- 
parameter Galster fit given in [78]. The Kaskulov and 
Grabmayr model has the same physical basis as that of 
Priedrich and Walcher, but is much more constrained; 
nevertheless, it fits the Geu data quite well. This result 
suggests that the radius of the 3q nucleon core is 

(r-')3f = ^ = (0-68±0.01)fm. (29) 



3. Quark models 

The predictions of several recent relativistic con- 
stituent quark models are compared with the data in Fig. 
29. All employ a linear confining potential. The solid 
curves show calculations of the Pavia-Graz collaboration 
[110] that used the point-form spectator approximation 
(PFSA) for pointlikc constituent quarks and a Goldstone 
boson exchange interaction fitted to spectroscopic data. 
No additional parameters were adjusted to fit the form 
factors. The data for Gep/Gd are reproduced very well 
and the data for magnetic form factors are also described 



relatively well for < 1 (GeV/c)^, but the calculated 
value of Gmp/^pGd decreases too rapidly for larger Q^. 
The prediction for Geu/Gd lies well below the data for 

> 1 (GeV/c)^. The long-dashed curves show calcu- 
lations of Simula [111], based upon the model of Car- 
darelli et al. [112], that used the light- front approach 
and the one-gluon exchange (OGE) interaction. Here, 
constituent-quark form fac;tors were fitted to data for 

<1 (GeV/c)^ and the calculations were extrapolated 
to larger . This approach provides good fits up to 
about 4 (GeV/c)^. Finally, the dash-dotted curves show 
the results for a semirelativistic hypercentral constituent 
quark model (hCQM) [113] where the constituent-quark 
form factors, chosen as linear combinations of monopolc 
and dipole forms, were also fitted to recent data. Of the 
selected quark model calculations, their results clearly 
achieve the best overall agreement with the data. 

Finally, the most recent lattice QCD calculations of 
nucleon form factors were reported by the QCDSF col- 
laboration [114] using nonperturbatively improved Wil- 
son fermions in the quenched approximation. Unfortu- 
nately, straightforward chiral extrapolation [115] does 
not provide adequate agreement with data for < 
1.5 (GeV/c)^. Matevosyan et al. [116] proposed a 
model-dependent extrapolation procedure based upon 
the LFCBM. This extrapolation is quite severe because 
the lattice calculations remain limited to quark masses 
that correspond to > 0.5 GeV, lattice spacings with 
a > 0.05 fm, and volumes that might not fully contain 
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the pion cloud; therefore, comparison with data is prob- 
ably premature. 



VII. SUMMARY AND CONCLUSIONS 

We reported values for the neutron electric to magnetic 
form factor ratio, GEn/GMni deduced from measure- 
ments of the neutron's recoil polarization in quasielas- 
tic ^H(e, e'n)'^H kinematics at three acceptance-averaged 
Q2 values of 0.45, 1.13, and 1.45 (GeV/c)^. In the 
one-photon exchange approximation for elastic scattering 
from a free neutron, the polarization vector of the recoil 
neutron is confined to the scattering plane and consists of 
a longitudinal component, Pf^^ oc and a transverse 

component, P^^'^ oc GsnGMn- The use of a deuteron 
target to access the neutron form factor ratio via the 
quasielastic ^H(e, e'n)^H reaction has the advantage, as 

established by Arenhovel et al. [18, 19], that both P^^^ 
and P]''^ are relatively insensitive to final-state interac- 
tions (FSI), meson-exchange currents (MEG), isobar con- 
figurations (IC), and to theoretical models of deuteron 
structure. 

A high-luminosity neutron polarimeter designed specif- 
ically for our experiment, Jefferson Laboratory EOS- 
OSS, was used to measure neutron polarization-dependent 
scattering asymmetries proportional to the projection of 
the polarization vector on the transverse axis. A dipole 
magnet located upstream of the polarimeter was used to 
precess the neutron polarization vector in the transverse- 
longitudinal plane, thereby permitting access to the ra- 
tio Pj'^^/Pj'''' oc GEn/GMn- Values for the scattering 
asymmetries were extracted from neutron time-of-flight 
measurements in our polarimeter via the cross ratio tech- 
nique. The merit of the cross ratio technique is that the 
scattering asymmetries are independent of the luminosi- 
ties for the two electron beam helicity states, and inde- 
pendent of the efficiencies and acceptances of the top and 
bottom halves of the polarimeter. Systematic uncertain- 
ties in our results are minimal as the analyzing power of 
the polarimeter and the polarization of the electron beam 
cancel in the form factor ratio. Further, other sources of 
uncertainty, such as radiative corrections and neutron de- 
polarization by lead shielding, are small as they nearly 
cancel in the ratio. 

To account for the finite experimental acceptance and 
nuclear physics effects (i.e., FSI, MEC, and IC), we used 
two independent simulation programs to average theo- 
retical ^H(e, e'n)^II recoil polarization calculations com- 
puted according to the model of Arenhovel et al. [18, 19] 
over the acceptance. The results from these two sim- 
ulation programs agreed to better than 1% at our two 
lower points and 2% at our highest point. Fur- 
ther, by averaging two difi'erent sets of theoretical calcu- 
lations assuming different parametrizations for Geu, our 
acceptance-averaged and nuclear physics corrected values 
for Geu were found to be insensitive to the choice of the 



dependence of Geu- 

Our results for Geu and data on the nucleon form 
factors were compared with selected theoretical model 
calculations. All of the model calculations based upon 
vector meson dominance and those emphasizing the pion 
cloud presented here provide qualitative agreement with 
some of the four nucleon form factors, but no model 
achieves simultaneous agreement with all four form fac- 
tors. The predictions of several recent relativistic quark 
models also achieve qualitative agreement with the data, 
with the most successful models utilizing form factors 
for the constituents; the results from a chosen model as- 
suming pointlike constituents are not as successful. Al- 
though a comparison between data and the results of lat- 
tice QCD calculations is probably premature, the recent 
precise data obtained from experiments employing polar- 
ization degrees of freedom will no doubt serve as a future 
challenging test of QCD as formulated on the lattice. 

In conclusion, our results at = 1.13 and 1.45 
(GeV/c)'^ are the first direct measurements of Geu using 
polarization degrees of freedom in the > 1 (GeV/c)^ 
region and are the most precise determinations of Geu 
over all ranges of . The achievement of relative statis- 
tical uncertainties in the form factor ratio GeuIGmu of 
8.4% and 9.5%, respectively, at these two points, to- 
gether with relative systematic uncertainties on the level 
of 2%, was a triumph for our high figure-of- merit and 
high luminosity neutron polarimeter. 
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APPENDIX A: FORMALISM FOR THE 
QUASIELASTIC ^H(e,e'n)^H REACTION 

Our notation for the kinematics and nucleon recoil po- 
larization for the quasielastic ^H(e, e'n)^H reaction fol- 
lows that of Arenhovel et al. (e.g., [117]). For ease of 
notation, all kinematic quantities in the center-of-mass 
(cm.) frame of the recoiling neutron-proton (n-p) system 
will carry a superscript '^•™'; however, kinematic quanti- 
ties referred to the laboratory frame will not be adorned 
with a superscript 
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FIG. 30: Schematic diagram of the kinematics for the electrodisintegration of the deuteron in the one-photon exchange approx- 
imation as viewed from the laboratory frame and the recoiUng n-p cm. frame. 



Kinematic notation 



where 



A schematic diagram of the kinematics for the elec- 
trodisintegration of the deuteron in the one-photon ex- 
change approximation is shown in Fig. 30. Our no- 
tation for the electron kinematics is as usual, and we 
assume the electron scatters from a deuteron with ini- 
tial four-momentum (md,0). Following the break-up of 
the deuteron, the proton and neutron exit with three- 
momenta Pp and p„, respectively. As is customary, we 
use 6pq (Onq) to dcuotc the polar angle between Pp (p„) 
and q in the laboratory frame, and a reaction plane is 
defined by any two of q, Pp, and p„. As is shown in 
Fig. 30, the reaction plane is tilted at a dihedral angle <j) 
with respect to the scattering plane. It should be noted 
that in the n-p cm. frame, this dihedral angle, ^!>^'™', is, 
obviously, just equal to (f). 

The n-p cm. frame is reached via a boost along q. In 
the laboratory frame, the n-p final state has an invariant 

q^, where the relative n-p 



mass. 



W„p, ofWnp = 



^2 

np 



energy in the laboratory frame, £„p, is Enp = u> + ma- 
With these definitions, it is clear that the Lorentz factor 
for the boost from the laboratory frame to the n-p cm. 
frame is 



En 



np y/{uJ + TOd)2 - q2 



(Al) 



We denote the polar angle between the relative n-p mo- 
tion in the cm. frame, pj;;™' = ^(Pp™' — Pn™') = 
Pp™' (assuming equal nueleon masses), and q'=™- as 
Qnp^'- As can be shown easily, this angle can be writ- 
ten sole ly in term s of the laboratory frame observables 
En = \/pl + ml, \pn\, 9nq, and u) as 



cose-p-- = 



|Pn| cos6>„g - \q\E„/E„ 
VA + B 



(A2) 



np > 



h\En 



B = |p„| COsOnq 



E.^ 



np 



(A3a) 
(A3b) 



Clearly, G^p"^- 



0° (180°) corresponds to perfect 
quasifree emission of the proton (neutron); however, it 
should be noted that there is vanishing phase space for 
perfect quasifree emission. 



2. Recoil polarization 

The five-fold differential coincidence cross section for 

the electrodisintegration of the deuteron in polarized- 
electron, unpolarized-deuteron scattering is of the simple 
form [19] 



cr(/i,0,0) = 



dV 



diJe'dfie'dfi^-^- 



= ao (1 + hP^A^) , (A4) 



as the electron asymmetry, A^, is the only polarized con- 
tribution to the cross section. As usual, (Tq denotes the 
unpolarized cross section. The above expression for the 
cross section can also be written in terms of structure 
functions as [19] 

cj{h, 0, 0) = C{pLfL + PtIt + PltUt cos <^^p"- 

I „ r o J,c.m. 



+ ptt/tt cos 24)';;™' 

+ hPep'LTf'LT^^^4>np ), 



(A5) 



where the structure functions are evaluated in the n- 
p cm. frame, the pi are elements of the virtual photon 
density matrix and functions of kinematics, and C is a 
function of kinematics. It should be noted that the above 
expression for the cross section is differential in E^i , fig' , 
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and O^-^-. The Jacobian, J = 9f7^-^7an„, which trans- 
forms fijj™' is given by [117] 



J = - 



(inln 



^ \ /?c.m.-,c.m. 
' \t-^n In 



cos(7r — O! 



c.m. 
np 



(A6) 



Here, 7 is as given in Eq. (Al), 7^™' is the Lorentz factor 
for the boost which takes the neutron from its rest frame 
to the n-p CM frame, 



In 



2m J 



(A7) 



and 7„ is the Lorentz factor for the boost which takes 
the neutron from its rest frame to the laboratory frame, 



In = ll'n"^- [1 + m"^- C0s(7r - e-"-)] 



(AS) 



/3, and /?„ are the vefocities associated with 7, 

7^'™', and 7„, respectively. 

The nucleon recoil polarization in the n-p c.m. frame, 
P'^-™-, is of the form [19] 



dV 



di;e'dOe'dO^-^- 



= ^o[(p(o))-- + /,p^(pW)' 



(A9) 

where P^"-* and P*^'*' denote, respectively, the helicity- 
independent and hclicity-dependent recoil polarization. 
Written in terms of g''"'^ structure functions, the helicity- 
independent polarization components are 



C 



(p(O)^c... 



C 

— [plqI + Pt9t + PltqIt cos <j)\ 

CTO 



c.m. 
np 



+ /9tt5?t COS 2Cp™-)' 



np 

(AlOa) 
(AlOb) 



C_ 



^Plt91t sin + PttQtt sin 20^, '"- 



(AlOc) 

and the helicity-dependent polarization components are 

-,(/()xc.m. _ C 



(Plt5lt cos 4>l^- + p'tq't) , 



C 



,c.m. 
rip ' 



(p(/.))C.m. ^ ^ (pi^^iVcOSc/,^-- 



(Alia) 

(Allb) 

(Allc) 



The boost from the laboratory frame to the n-p c.m. 
frame is along q, which is not, in general, parallel to 
either nucleon's momentum vector; therefore, the recoil 
polarization components in the laboratory frame arc re- 
lated to the recoil polarization components in the n-p 



c.m. frame via a relativistic Wigncr spin rotation. As the 
nucleons' momenta span the t-i plane, the n-component 
is unchanged, while the t- and .^-components mix accord- 
ing to 



3C.in. 

■ 3 ' 



(A12) 



where i,j G {t, n, £}, TZij (0^) denotes a matrix element 
of the Wigner rotation matrix, 



cos0^ sin6'Jf 

1 
-sin6l!f coseW 



(A13) 



and 6^ , the Wigner rotation angle for the neutron, is 
expressed in terms of kinematics as [117, 118] 



= sin-i 



1+7 



■sin(e^----e„) 



(A14) 



Here, ^5j™'(= tt— Qnp^ ) and On denote, respectively, the 
polar angle of the neutron's momentum vector relative 
to q in the n-p c.m. frame and the laboratory frame. For 
non-relativistic boosts (i.e., 7, 7^'™', and 7„ all ~ 1), it 
is clear that we recover the non-relativistic result, 0^ — > 
6*^™- — On- Also, it is obvious that for perfect quasifrcc 
emission (i.e., O^™' = or tt), the recoil polarization 
components in the n-p c.m. frame are identical to those 
in the laboratory frame. 



APPENDIX B: SENSITIVITY TO NUCLEAR 
PHYSICS EFFECTS AND DEUTERON 
STRUCTURE 

To demonstrate the sensitivity of P^^^ to the value 
of Geu and the inscnsitivity of p/''-' and pj^'^^ to FSI, 
MEC, IC, and the choice of the NN potential, we 
present several examples of ^H(e, e'n)^H recoil polar- 
ization calculations performed within the PWBA and 
FSI+MEC+IC+RC models of Arenhovel et al. [18, 19, 
119] in Figs. 31, 32, and 33. We have (arbitrarily) cho- 
sen to show examples of these calculations for the cen- 
tral kinematics of our Q"^ = 1.136 (GeV/c)^ point (i.e., 
Ee = 2.326 GeV, E^' = 1.718 GeV, 6'e' = 30.93°). 

First, FSI+MEC+IC+RC calculations of p/'*^ arc 
shown in Fig. 31 as a function of Q^p^' for three values of 
Gsn scaled by the Galster parametrization: 0.5, 1.0, and 
1.5. A strong (nearly linear) sensitivity of Pj;^^ to the 
value of Gsn is seen at and near quasifree emission. Sec- 
ond, the inscnsitivity of P^^^ and P^^^ to FSI, MEC, and 
IC for quasifree emission is shown in Fig. 32 where little 
difference between the PWBA and FSI+MEC+IC+RC 
calculations is observed at and near quasifree emission. 
Finally, we compare FSI+MEC+IC+RC calculations of 
P/'*^ and P/'*^ for the Argonnc VIS [120], Bonn [74], 
Nijmcgcn [121], and Paris [122] NN potentials in Fig. 
33. Again, at and near quasifree emission, there is little 
model dependence. 
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FIG. 31: (Color online) Sensitivity of FSI+MEC+IC+RC cal- 
culations of Pf'^^ to the value of Geu for the central kinematics 
of our — 1.136 (GeV/c)'^ point. The results shown are for 
(t>np^' ~ 0° and the Bonn potential. 
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FIG. 32: (Color online) Comparison of PWBA and 
FSI+MEC+IC+RC calculations of P/''' (top panel) and P^''"' 
(bottom panel) for the central kinematics of our Q2 = 1.136 
(GeV/c)^ point. The results shown are for (?f>5i™' = 0° and the 
Bonn potential. 



of 



0.3 


A 


Argonne VI 8 


0.2 


\\ 


Bonn 

Nijmegen 


0.1 


\ 


Paris 



0.1 






J \ 
\ 


/ — ^%\ -. 

<^ \J- 


0.2 






0.3 







20 40 60 80 100 120 140 160 180 



0.6 - 




20 40 60 80 100 120 140 160 180 



[deg] 

np 

FIG. 33: (Color online) Comparison of FSI+MEC+IC+RC 
calculations of P/''' (top panel) and P^'''' (bottom panel) for 
the Argonne V18, Bonn, Nijmegen, and Paris potentials. The 
results shown are for the central kinematics of our = 1.136 
(GeV/c)2 point and c^^'^' = 0°. 
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